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1.0 INTRODUCTION 


V SAERO is a computer program for calculating the non-linear 
aerodynamic characteristics of arbitrary configurations in sub- 
sonic flow. Non-linear effects of Vortex Separation and Vortex-' 
Surface interaction are treated in an iterative wake-shape calcu- 
lation procedure, while the effects of viscosity are treated in 
an iterative loop coupling potential flow and integral boundary 
layer calculations. The program is under continued development 
(Refs. 1 through S) in a number of directions outlined in Figure 
1. The present document is an interim User's Manual for the 
basic non-linear capability identified in Figure 1, The program 
scope covered by this document is given in Section 2.0. 

The basis of the computer program is a surface singularity 
panel method using quadrilateral panels on which doublet and 
source singularities are distributed in a piecewise constant 
form. The panel source values are directly determined by the 
external Neumann boundary condition controlling the normal 
component of the local resultant flow: the doublet values are 
solved after imposing the internal Dirichlet boundary condition 
of zero perturbation potential at the centers (underside) of all 
the panels simultaneously. Surface perturbation velocities arc 
obtained from the gradient of the doublet solution, while field 
velocities are obtained by direct summation of all singularity 
panel contributions. An outline of the mathematical formulation 
for the basic method is given in Section 3.0. 

Details of the configuration modelling are given in Section 
4.0, while Section 5.0 shows how the off-body velocity scan 
volumes are formed. Section 6.0 describes the way the program is 
run on the C.D.C. Cyber 176. 

The program input is described in Section 7.0 in three ways: 
first, the input variables are listed card by card; secondly, the 
function of each variable is described in detail on a card by 
card basis, and finally, a flow chart of the input is given 
showing in a more visual way the sequence and relationship of 
cards according to the various input options. This latter form 
is especially useful when some familiarity with the input 
variables has been gained. 

The program output is described in Section 8.0, while the 
input and output for an example case of a wing-body configuration 
is given in the Appendix. 

The program is written in standard FORTRAN IV and has been 
developed on the C.D.C. Cyber 176 computer. Minor changes to the 
code allow it to run on IBM and Cray computers. There are two 
versions of the basic code; the first allows up to 1,000 panels 
(i.e., unknowns) and the second allows up to 3,000 panels. Where 
limiting values for variables are quoted the 1,000 panel numbers 
comes first followed by a slash then the 3,000 panel number. For 
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example, the core requirement (octal) for the program is 
234K>203K in SCI! and 0-^3 2 IK in LCU. 

The CPU time required by the program is plotted in Figure 2 
against the number of panels (unknowns). The CP seconds are 
converted to CDC 7600 values. The run times are for cases with 
symmetry about the x - z plane; i . e . , the actual number of panels 
is double the value shown. The approximate effects of wake 
complexity and iteration cycles are indicated in Figure 2. 
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Figure 2. CPU Time (CDC 7600 C.P. secs.) 


2.0 PROGRAM SCOPE 


This section describes the general capability of the program 
covered by the present user manual. Other options are in the 
code--in fact, several of these appear in the input lists here; 
however, these are still in the process of being checked out or 
are being developed further. 

The program was generated primarily for analyzing high-lift 
configurations; however, a number of additional features now 
allow the program to be applied to complete aircraft 
configurations, including wings, bodies, tailplanes, fins, slats, 
slotted flaps, powered nacelles, etc. The non-linear wake 
relaxation routines for the high-lift calculations are applicable 
also to the full configuration cases. 

The analyses are normally performed in 'symmetrical' flight; 
i.e,, only one half of the aircraft geometry is described. 
Asymmetric conditions can be treated but the eoaglete geometry 
must be described at this time (i.e., even if the geometry has 
ay suae try ) . 

The configuration can be analysed in ground effect; the 
ground plane is actually the horizontal x.-y plane of the global 
coordinate system. The components of the configuration can be 
moved and rotated within the program to set up the required 
conditions. Non-zero normal velocities can be specified on user- 
selected sets of panels to represent inflow or outflow; e.g., for 
powered nacelles. Vortex sheet wakes, representing the shear 
between outflow (jet) velocity and local external velocity, are 
then attached to the surface panels to enclose each jet. 

Although the program deals primarily with closed surfaces 
(i.e., t.hi.c.k wings and bodies), provision is made to treat open 
surfaces also. In such cases only the Neumann boundary condition 
is satisfied and the singularity model for that part of the 
surface is identical to that used in the quadrilateral vortex 
method (Ref. 6). Open and closed surfaces can be mixed in a 
configuration. 

Surface streamline calculations can be requested. Each 
calculation starts in a panel and proceeds upstream and 
downstream until either a stagnation region is entered or until 
the calculation reaches an edge or a cut in panel neighbor 
relationships (see 3.2). The user supplies a set of 'starting' 
panels to adequately cover the region(s) of interest. 

Two integral boundary layer calculation options are provided 
in a viscous-^potential flow iteration coupling. The user 
specifies the number of iteration cycles. 

At the end of the calculation the user can request off-body 
velocity surveys and off-body streamlines. 



The program generates a plot file of geometric and 
aerodynamic data® The user must save this in the JCL (see 6.2) 
for plotting at a later stage® The program can also generate 
restart files for the purpose of continuing with the calculation 
at a later stage (see 6*3)* Provision is made to change some of 
the input data on a restart run. 
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3.0 MATHEMATICAL MODEL OVERVIEW 

3 • 1 Formulatio n 

The general arrangement of a configuration is shown in 
Figure 3 and is described in terms of a body-fixed GLOBAL 
COORDINATE SYSTEM referred to as the G.C.S. The configuration is 
immersed in a uniform onset flow, ^ with velocity potential ^ . 
Potential flow is assumed to exist both inside and outside the 
boundaries of the configuration. In the external flow field 
(i.e., the one of interest) the total velocity potential, 4> , is 
the sum of the onset flow potential, ^ , and the perturbation 
potential, $ . Similarly, the total velocity potential insid e the 
configuration is For the present description, the wake 
surfaces are assumed to have vanishing thickness with zero 
entrainment . 

After applying Green's Theorem to the inner and outer 
regions and combining the resulting expressions, the velocity 
potential at a point P on the inside surface can be written 



47T$ d = f f - V" * V (r) dS " 27r($ " V p 

S-P 

+ // (^u-\)n- v(l)dW 
W 


+ 



(v$. - v$) dS + 4^ P 


(l) 


where r is the length of the vector from the surface element to 
the point P, and S-P signifies that the point P is excluded from 
the surface integral. $ u - $ i is the local jump in potential 
across the wake surface, W. Equation (1) gives the total poten- 
tial at the interior point, P, as the sum of perturbation poten- 
tials due to a normal doublet distribution of strength 
( $ - $ 4 ) on S and ( $ tt - on respectively, and a source 

distribution of strength, - V$) ©n S. The potential for 

the uniform onset flow, (j)^ , is also included. 

In principle, an infinite number of combinations of doublet 
and source distribution will give the same external flowfield, 
but different internal flowfields. To render a unique 
combination of singularities either one of the singularity 
distributions must be specified (e.g., o * 0 in the doublet-only 
formulation) or, as in the present case, the internal flow must 
be specified. There are several reasonable options for the 
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internal potential flow. On the of most convenient options is to 
set = 4 a,- This internal Dirichlet boundary condition gives 
zero perturbation to the onset flow inside the configuration. 
Eqn . (1) becomes: 


0 = 



<J>n • V 


(*U ’ 


(1) dS - 2n* p 
\)n ■ V (1) dW 
(V^ - VS>)dS 


( 2 ) 


where ii» the perturbation potential in the flow field, has been 
substituted for $ - $ . 

OO 

The first two terms in Eqn. (2) give the perturbation poten- 
tial due to a distribution of normal doublets of strength, 4, on 
the configuration surface, S. Similarly, the third term rep- 
resents a doublet distribution of strength, $ u - $ L* on wa ^ e 

and the fourth term represents a source distribution of strength 
n • ( VcJ)^ - V3>) on the configuration surface. 

For the problem of analyzing the flow about a given 
configuration geometry, the doublet distribution, , is the 
unknown while the source distribution is determined directly by 
the external Neumann boundary condition: 

n * V$ = -V N 

where Vn * s the local normal velocity component of the external 
flow relative to the fixed surface. (This component is positive 
in the direction of n.). 

For solid boundary conditions, i s ze ro; for more general 

cases, however, Vn can l iave several parts representing various 
special effects (Refs. 1 through 5). Two of these effects are 
considered here: 

(i) boundary layer displacement effect using the 
transpiration technique; and 

(ii) i nf 1 ow^out f 1 ow for engine i n 1 e t^e xhau s t representation. 
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Thu s * 


V N = ys< u e<5*) + VnoEU 

where the first term represents the rate of 'growth 1 of the 
boundary layer and the second term, ^NORM* is positive for out- 
flow and negative for inflow. The boundary layer term is usually 
zero at the start of the calculation and is then updated during 
each v i scous^po t ent ia 1 iteration cycle. 

Thus the source strengths can be evaluated in Eqn. (2) 
before each doublet solution: 


(V<fc - V4>) = 4 TO = 


3s 


(U e 6*) 


+ V 


NORM 


- n 


1 


The above internal Dirich let boundary condition is 
applicable for closed surfaces only. If there are parts of the 
surface that are extremely thin (e.g., thickne ss-^chord ratio < 
5 ( .») , or are wing-like and are well away from the area of 
interest, then these can be represented by open sheets using the 
Neumann boundary condition only. The treatment of these open 
surfaces is identical to that in the quadrilateral vortex method 
(Ref. 6). The thin and thick options may be mixed within the 
matrix of influence coefficients and so a configuration can 
include patches of both types. 


3.2 ur £ 


Figure 4 shows a flow chart of the numerical procedure for 
the complete solution. For this purpose the configuration sur- 
face and wake are represented by a number of flat quadrilateral 
panels (Figure 3). (The way the panels are formed will be 
treated in Section 4.) The surface integrals are performed in a 
piecewise manner with the assumption of uniform singularity 
distributions over each panel. Eqn. (2) is then satisfied simul- 
taneously at a control point on each surface panel, i • e • , Eqn. 
(2) becomes 


£ (p K C JK> - 2 T] + 
K=1 



J=l, NS 

(3) 


where N s , N w are the number of surface and wake panels, 
respectively. The panel doublet value (unknown) is d j' 4 jt . The 
quantities Bj j' , Cjjt are the perturbation velocity influence 
coefficients for the constant source and doublet distributions, 
respectively, on panel K acting on the control point on panel J. 
(The control point is the average point of the panel's four 
corners.) These coefficients include contributions from 'image' 
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Figure 4. VSAERO Program Outline. 






for cases of symme try and/or giosud effect. (For control points 
on open surfaces the Neumann boundary condition equation has a 
similar form to Eqn. (3), but the influence coefficients are for 
the normal component of induced velocity.) 

Equation (3) is solved by a direct method when N # £ 320^630 
and by a blocked Gauss-Seidel iterative procedure for N $ value s 
above this. 

The surface gradient of p is evaluated on each panel by 
differentiating a two-way parabolic curve fit through the doublet 
value s on the panel and its four immediate neighbors. At certain 
lines on the surface where there is a jump in local conditions, 
e.g., wake separation lines, edges of inf low-'outf low regions, 
etc., the program 'steps back 1 and uses information from the 
'neighbor of a neighbor' to avoid taking a gradient across such a 
line. The doublet gradient provides the local tangential 
velocity perturbation which is combined with the local tangential 
component of the onset flow to give the total local velocity. 
This is used to evaluate the pressure coefficient at each panel's 
center. The integrated force and moment for the configuration 
and its parts are obtained by summing the pressure force and 
moment contributions from each panel. 

Then a solution has been formed, a wake shape iteration loop 
can be executed, In this loop the wake is repositioned so that 
the streamwise edges of wake panels are aligned with the local 
calculated flow directions. The wake panels and their influence 
coefficients are then reformed prior to a new solution. Figure 4. 

When the wake shape iteration is completed the 
viscouso'potential iteration loops can be entered. At this time 
there are two integral boundary layer options. The first (B.L.l 
in Figure 4) is an infinite swept wing stripwise calculation 
intended primarily for wing-type surfaces. This method is des- 
cribed in Reference 7. It uses calculated pressure and velocity 
distributions from streamwise strips on the wing. The second 
option, B.L.2, is a two-dimensional boundary layer calculation 
along calculated streamlines; local curvature effects are treated 
as if the surface were part of an axisym-metric body (Ref. 8). 
Use of the B.L.2 option requires surface streamline calculations 
to be performed and is applied mainly on body-type surfaces. 

Both boundary layer routines return transpiration source 
values back to the potential flow code. Since the source values 
only affect the right-hand side of the equations, a new solution 
can be readily obtained (Figure 4). The boundary layer routines 
also provide the locations of separation lines; if there is a 
large amount of separation or if the separation zone ha s changed 
markedly from the previous iteration solution, then a new wake 
must be formed before the next solution. This part of the pro- 
cedure is still being checked out. 



After the last iteration* off-body velocity surveys and off- 
body streamlines may be executed. These use a general velocity 
calculation scheme somewhat similar to Eqn* (3) except velocity 
influence coefficients are used instead of velocity potential 
coefficients* Special near- field techniques are activated for 
velocity points that are close to the body surface or wake* 
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4 . 0 CONFIGURATION MODELLING 

4 . 1 O vervie w 



This section describes the way the configuration is 
subdivided to form the set of panels used in the numerical 
solution. The operation has been designed with a view to both 
user convenience and to internal efficiency in the program code. 
First, the panels are organized in a multiple level hierarchy to 
facilitate treatment of complicated configurations and also to 
simplify relocating separate pieces of a configuration (e.g., 
flaps): inside the code the organization simplifies the problem 

of identifying a panel's immediate neighbors. Secondly, options 
for automatic paneling are provided so that it is not necessary 
for the user to provide the actual panel corner points. (He 
still £a_u do this if he wants to.) T h ese options also provide a 
powerful capability for changing, in a later run, the emphasis of 
high panel density from one part of the configuration to another. 
Such a repaneling is achieved with only minor changes to the 
input file, the input geometry can remain essentially the same. 


4 . 2 Hi erarch y 

The description of the mathematical model in 3.1 has already 
shown the initial breakdown of a CONFIGURATION into two major 
parts; viz., the (solid) SURFACE and the (flexible) WAKE. In 
fact, on a complicated configuration, there may be several 
separate surfaces and several wakes. Each wake is broken down 
directly into columns of WAKE PANELS while the surface is 
subdivided into several parts; i.e., COMPONENTS, ASSEMBLIES, 
PATCHES and PANELS, Figure 5. (Panels are further subdivided 
into SUBPANELS, but this is an internal automatic feature 
activated for special treatment of near-field velocity 
calculations . ) 

The smallest 'stand-alone' part of a configuration surface 
is a PATCH. A patch contains a rectangular array of PANELS. 
Patches are grouped together under two headings (Figure 5): 
COMPONENTS and ASSEMBLIES. These are at the same level in the 
hierarchy because they each contain one or more patches, in fact, 
they overlap. The main difference between the two groupings is 
that the patches within an assembly are regarded as detached from 
patches on other assemblies while the patches within a component 
may or may not be contiguous with patches of other components. 

The nature, purpose and treatment of each part of the 
configuration are described in the following subsections; in the 
meantime. Figure 6 shows an example of the breakdown of a high- 
lift wing configuration into assemblies, components and patches. 
The wing has a leading -edge slat and part-span slotted flap. 
Thus there are three separate pieces of surface which are 
identified as separate assemblies. For this example, the wing 
and slat have been placed under component 1 and the flap is 
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j^SUBPANELSj [SUBPANELSj 


Figure 5. Configuration Hierarchy. 
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Of 

ASSEMBLY 1 ASSEMBLY 2 



ASSEMBLY 3 
(FLAP) 


Figure 6. Examples of Assemblies, Components and 
Patches on a High-Lift Configuration. 



called Component 2 — this grouping is entirely up to the user. 
The grouping shown allows the flap to be easily repositioned to a 
new setting later on (see the component subsection). Typical 
patches (A, D , C , D ) are indicated in Figure 6. A s will be 
de scribed below in 4. 5. 7. 2, automatic procedures can be called on 
to form the t ip-edge patches closing the wing, flap and slat. (A 
small closing patch can be left off when a detailed solution is 
not required in that area.) 

4 . 3 As s emb ly 

An assembly is a completely separate piece of the surface. 
For example) the entire wing, fuselage, fin, tailplane and 
nacelle of an aircraft configuration would form one assembly, 
while slotted flaps, slats, etc., would be separate assemblies. 
(Small attachment structures would not normally be represented in 
the global treatment of a full configuration.) 

The primary function of the assembly status is to prevent 
panel neighbor relationships being formed across small slots or 
gaps. Each panel needs to 'know' its four immediate neighbors 
for the purpose of evaluating the doublet gradient (3.2). Within 
a patch the neighbor relationships are straightforward, but 
across patch edges the panel neighbors are identified by an 
automatic procedure. During its search for possible neighbor 
candidates, this procedure only considers panels within the same 
assembly. 

A secondary function of an assembly is that it has its own 
printout of integrated force and moment data. 

The assembly number for each patch is identified on the 
patch card (CARD 10) using the variable KLASS. 


4 . 4 C omp o ne nt 

The set of patches forming a component is completely in the 
hands of the user and the patches need not necessarily be 
contiguous* The main function of a component is to provide the 
user with a separate printout of integrated force and moment data 
for his selected parts of the configuration. The component 
number for a patch is set by an internal counter. This counter 
is incremented by one when the user places a value of 4 for NODES 
on a patch's terminating section card (CARD 11). (Note that a 
component card (CARD 9) for the new component must then precede 
the next patch card (CARD 10)). 

It should be remembered that the order of patches in the 
input is dictated more by surface proximity (see 4.5.4) rather 
than the user's choice of patch sets for a component. For this 
reason, provision is made to override the current component 
counter using the variable ROMP on the patch card (CARD 10). 
This does not change the counter value, it merely identifies the 
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present patch with another component; this facility must not be 
used oi% the f j.r s t patch of a new component. (Why change the 
counter otherwise?) 

Each component can be defined in its own CO li POM ENT 
COORDINATE SYSTEM (referred to as the C.C.S.). This facility 
allows all the patches within a component to be scaled, moved and 
rotated (in that order) en bloc with one instruction on the 
component card (CARD 9). 

The component card contains the appropriate transformation 
information which converts the geometry from the C.C.S. to the 
global coordinate system (referred to as G.C.S.). This 
information includes (i) the translation vector, (CTX, CTY, CTZ), 
which is simply the origin of the C.C.S. located in the C.C.S. ; 
(ii) the scaling factor; and (iii) the rotation angle, 6, about a 
hinge line vector, h, Figure 7. Provision is made for the user 
to specify two points on a general hinge line vector (in the 
C.C.S.), otherwise the y-axis in the C.C.S. is used. Both the 
scaling and the rotation are applied in the C.C.S. prior to the 
translation. This component transformation is performed at the 
end of the geometry input routine; i.e., after the basic geometry 
of the complete configuration is assembled. 


4 . 5 Pa tch 

4.5.1 Sh si££ 

The developed (i.e., 'opened out') shape of a patch should 
be roughly four-sided to keep panel shapes and distributions 
reasonably regular. This does not exclude the presence of kinks 
in any or all of the patch sides, but kink angles should not be 
large (the upper limit has not been established, but for the time 
being 6C° should be regarded as a large kink angle). One side or 
two opposite sides of a patch may be reduced to zero length 
provided the overall patch shape is reasonably regular. Figure 8 
gives some basic guidelines for acceptable patch shapes. 

Awkwardly shaped surfaces may be represented by several 
patches, whereas simple shapes may require just one patch. In 
Figure 6, for example, parts A, B , C and D are typical patches on 
assembly 2 . The main surfaces of assemblies 1 and 3 may be 
formed by single patches; additional patches may be used to cover 
the open ends (shaded). Patches that are wrapped around until 
two opposite sides meet (e.g., patches A, B and C in Figure 6) 
are referred to as FOLDED PATCHES. It is recommended that folded 
patches (with sides 2 and 4 meeting) be used to represent all 
wing-type surfaces. 
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Figure 7 . Component Transformation. 
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AVOID LARGE KINK ANGLES — USE THREE PATCHES HERE 
Figure 8. Preliminary Guidelines for Patch Shapes. 
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4.5.2 Convention 


Ve have seen above that a patch is basically a four-sided 
shape when developed or 'opened out'. 1 1 must always be regarded 
as such even if one of the sides or two opposite sides are made 
zero or if some of the sides have kinks. In the following 
discussions patches will often be regarded as rectangular --this 
is purely a convenience for discussion relationships and is not a 
shape restriction. Our view of the patch will always be from the 
outside ; i.e., looking onto the wetted surface from a point in 
the flow field. 

For convenience the terms 'chordwise' and ’spanwise' are 
used to describe the directions of the panel columns and rows, 
respectively. Figure 9. These directions are analogous to the 
conventional wing layout, but, in the patch context, these 
directions are not restricted to the x and y directions, 
respectively. For example, on a patch representing the wing tip 
it is convenient to have the columns of panels vertical and the 
rows of panels to be parallel to the wing chord; in this case, 
therefore, the 'chordwise' direction on the patch is actually 
vertical while its 'spanwise' direction is along the wing chord 
(see later in Figure 22). 

Patch geometry is defined using chordwise lines called 
SECTIONS. (These are described below.) A set of sections 
distributed spanwise across a patch defines the patch surface. 
The convention adopted here is that points defining a section 
shape proceed from top to bottom. Figure 10. (in the case where 
a patch represents the main surface of a wing, this convention 
causes the points defining each section to proceed from the 
trail ing-edge lower surface and finish at the tr a i 1 ing-e dg e upper 
surface.) In our view of the patch, the order of the sections 
proceeds in the positive spanwise direction. Figure 10. 

For the purpose of automatically connecting panels from one 
patch to another, it is important to identify patch sides. The 
convention adoped here is that the first and last sections 
defining a patch correspond to sides 1 and 3, respectively. With 
this convention, the order of the sides is a ntic lockwis e . Figure 
9. The order of the corner points follows the same sequence as 
the sides, starting with 1 at the top of side 1. All the panels 
within a patch take the same side and corner point convention as 
for the patch. 
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PANEL SEQUENCE 
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4.5.3 Parameters 

The main parameters associated with a patch are: 

PNAME, IDENT, KLASS, KOMP, NRQW, NCOL, I PAN, LPAN , NPAN 

PNAME is an optional text description for each patch while IDENT 
gives the patch type. At this time IDENT = 1 and 2 are used for 
patches on closed surfaces, 1 for a wing type and 2 for a body 
type. Both of these use the internal Dirichlet boundary 
condition. IDENT * 3 is for patches on open surfaces; on these 
the Neumann boundary condition is used. Different printouts are 
used for the analysis results depending on the IDENT value. 

KLASS and EOMP are the assembly number and component number 
associated with each patch and NROW, NCOL hold the patch's number 
of rows and columns of panels, respectively. I PAM and LPAN are 
the first and last panel subscripts (absolute) on the patch and 
NPAN is the number of panels on the patch. 

4.5.4 Input Ord er 

The order of input patches and their orientation is not very 
restrictive. The main requirement is that patches with panels 
'facing' each other across a small gap ( i . e . , small in relation 
to panel size) should be kept together in the input. The upper 
and lower surfaces of a wing in the trail ing -edge region should 
be part of the same patch to ensure that the upper and lower 
panels at the trailing edge occur in the same column of the 
patch; in this way they are sure to fall in the same matrix block 
in the Gauss-Seidel iterative solution procedure. 

Patches adjacent to each other and essentially in the same 
plane can be separated in the patch sequence--this is because a 
panel's potential influence coefficient at a point beyond its 
perimeter is small near the plane of the panel, 

4.5.5 Sec t ions 

4 . 5 . 5 , 1 S ect ion Co ordinat e Syst em 

Each section of a patch may be defined in its own local 
coordinate system, referred to as the SECTION COORDINATE SYSTEM, 
or S.C.S. The user provides the necessary information on the 
section card to transform from the S.C.S. into the C.C.S, This 
transformation is performed immediately a section's geometric 
description is complete. This transformation is separate from 
the component transformation in which the complete component is 
converted into the G.C.S. (at which stage the S.C.S. geometry is 
discarded). This double transformation — both levels of which are 
optional--offers useful flexibility when preparing the input 
data. One particular advantage is that the geometric 
relations hips--especially the rotations-- a re kept reasonably 
simple without sacrificing generality. 
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The information required to transform from the S.C.S. into 
the C . C . S . (see Figure 11) consists of: (i) the translation 

vector, ( STX, STY , STZ ) , which is the position vector of the 
S.C.S. origin expressed in the C.C.S. coordinates; (ii) a scaling 
factor which is applied in the S.C.S.; ( i i i ) the rotation angle 

( ALF , degrees) about the y-axis of t he S.C.S.; and ( i v ) the angle 
(THETA, degrees) in the C.C.S. x-y plane, between the projection 
of the S.C.S. y-axis and the C.C.S. y-axis. 

4 . 5 . 5 . 2 Ba s ic Points 

The contour line of each section is defined by a set of 
BASIC POINTS, (BX, BY, BZ). These points nay be used directly as 
panel corner points, i.e., MANUAL PANELING, in which case the 
user oust take care over the number of input points. 
Alternatively, an AUTOMATIC PANELLLNG ROUTINE, referred to as the 
A.P.R. nay be activated which interpolates through the basic 
points to form a new set of points corresponding to panel corner 
points. (Note, this is just a temporary set as the user nay opt 
to use the A.P.R. in the spanwise direction as well, in which 
case the sectio ns do not necessarily line up with panel edges.) 

No natter which paneli n g option has been selected, basic 
points should be reasonably dense in regions of high curvature, 
such as near the wing leading edge. 

Several options have beer, provided for defining the basic 
points and these, in combination with the two-stage 
transformation described above, provide great flexibility when 
preparing the input. The options may be exercised at the section 
level so the input form may be changed from section to section. 
The options available at this time are described below and are 
controlled by the value of INPUT. INPUT takes the value of 
IN MODE on the section CARD 11 for values of 1 through 4; these 
are illustrated in Figure 12, toget he r with instances for their 
use . 


INPUT values of 1, 2 and 3 are used when a section lies in 
one of the reference planes of the chosen S.C.S.; in these cases 
we have a constant coordinate, x, y, or z , respectively. With 
one coordinate fixed, we need input only two coordinates for each 
basic point, e.g., y and z when INPUT = 1. Provision is made to 
specify a third quantity to give a local adjustment to the 'con- 
stant' coordinate, e.g., when using INPUT = 2 we may specify x, z 
and 6y. Usually the 6-quantity is left blank (i.e., 0). The 
basic value for the constant coordinate is zero until the section 
points are transformed into the C.C.S., so the value of that 
coordinate in the C.C.S. must be provided in the transformation 
information on the section CARD 11, 

INPUT value of 4, which requires all the components of each 
basic point position vector, is used when defining a completely 
arbitrary se'etion shape. 
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VIEW IN DIRECTION OF S.C.S. Y- AXIS 


Figure 11. Section Transformation Into C.C.S. 
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An option to input basic points in polar coordinate format, 
see Figure 13, has been provided. This requires IN MODE = 1 2 on 
the section card. Other IN MODE opt ions have been provided for 
automatically generating basic points on certain shapes: these 

will be described below in 4. 5, 5. 3, 

Zero or negative INMODE values allow the present section's 
basic points to be copied over completely from any £r«v iou s.lx 
defined section. The section number is (-IN MODE) except when 
IN MODE = 0; the latter copies over the points from the section 
just completed. The section number specif if ed is the absolute 
number from the beginning of the input and includes other copied 
sections as well as sections which may have been generated 
automatically. If the section counting becomes complicated, 
alternative ways of copying are available as described under 
sspccial routines (4.5.7). The basic points are copied from the 
S.C.S, set ( i . e . , as originally specified) and are then 
transformed to the present C.C.S. according to the new section's 
transformation information, 

4 . 5 . 5 . 3 Cho rdwise Region s 

The basic points defining a section may be assembled in a 
number of CKORDV I SE REG IONS for the purpose of controlling the 
panel density and distribution on that section. In addition, the 
option on manual or automatic paneling is selected at the 
chordwise region level, allowing the user to switch from one to 
the other within each section wherever he' chooses. Chordwise 
regions are used only as an input co nv enience and are discarded 
in the program as soon as the surface paneling is complete, 

A chordwise region must end on a basic point called a NODE 
POINT, Figure 14. A NODE CARD, containing the chordwise region 
paneling information (see below), inserted after a basic p int in 
the input deck identifies that point as the end of a chordwise 
region. Mode points are usually placed at 'problem' areas where 
large velocity gradients are expected to occur, e.g., flap hinge 
line, leading edge, close-interference regions, but the user can 
place them wherever he wishes to change from one panel scheme to 
another or to ensure panel matching on patches facing each other 
across small gaps. Four types of node point are provided at this 
time and are described below. 

The information on a MODE CARD consists of just three 
integers. 

(i) NODEC identifies the node point and its type. 

(ii) NPC is the number of panels to be generated by the 
A.P.R. in the chordwise region just completed-~a 
zero value gives manual paneling. 


28 



Input: x, RAD, THET (degrees) 

Program computes: y « RAD * sin (THET) 

z = RAD * cos (THET) 


Figure 13. Polar Coordinate Input (INMODE=12) 





BASIC POINTS DEFINE 
SECTION SHAPE 



SEQUENCE OF 
BASIC POINTS 
ON EACH SECTION 


Figure 14. Chordwise Regions on a Section. 
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(iii) INTC controls the form of the distribution, in the 
automatic paneling mode and is inactive in the 
manual paneling mode. 

(The C on the end of each quantity distinguishes the 
chordwise from the corre sponding spanwise quantities, which end 
in S , 4.5.6) 

NQDEC values of 1 or 2 specify the end of a chordwise region 
with, respectively, continuous or discontinuous surface slope 
onto the next chordwise region. These values are, therefore, 
used only on regions ending in the interior of a section. The 

last point on a section is specified by NODEC = 3 and is the only 

node point that must always be specified even if manual paneling 
has been selected. Negative NODEC values are also permitted and 
initiate a special copying routine described in 4. 5. 7.1. 

Four panel spacing options are provided in the A.P.R. The 
action of INTC values of 0, 1 and 2 is illustrated in Figure 15, 
and is based on the cosine distribution giving increased panel 
density towards, respectively, the beginning and end, the 
beginning only, or the end only, of the region. Equal spacing 

throughout the region is provided by INTC = 3. Coupled with the 

flexibility offered by the choice of chordwise region location, 
these spacing options have proven adequate so far; however, other 
options could easily be added should the need arise later, e.g., 
one based on increments in integrated surface curvature, or on 
increments in doublet value from a preliminary two-dimensional 
solution for the section have been considered. 

Clearly, node cards provide the user with an extremely 
versatile paneling tool. With one card deck of basic points 
defining the configuration geometry, he can, from run to run, 
change the form of the paneling simply by changing two integer 
values on each node card. Not only that, he can also move node 
cards within the deck (but not the node cards at section ends) or 
remove some or add new ones from run to run. This allows the 
user to concentrate his paneling in areas of interest, leaving 
other areas more sparsely panelled. It thereby provides a very 
effective use of the limited number of panels available, yet, on 
a subsequent run a few small changes to the node cards allow the 
emphasis to be switched to another area without having to punch a 
new basic geometry card deck. 

There is just one important ground rule for the use of node 
cards; the Jto^al. number of panels (automatic and/or manual) on 
each section of a patch m us t b e the s ame. The total is, in fact, 
the number of panel rows, NHO W , for that patch. The program 
monitors the number of panels on each section and the calcula- 
tions are terminated with an error message should the user make a 
mistake. Provided this ground rule is satisfied, it is not 
necessary for the panel distribut ion to be the same from section 
to section --in other words, the number of chordwise regions and 
their node information can vary from section to section. 
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SURFACE DISTANCE FROM 



s l = SO - COS( 0|))/2 

WHERE 0j = (I-1)ir/N AND N IS THE NUMBER OF INTERVALS REQUIRED 

(<*) I NTC=0 



s, - S(1 -COS(0|)) 

WHERE 0| - (I - 1 )n/2N 

(b) IWTC=1 

Figure 15. Spacing Options 0, 



*j = S SIN(0j) 

WHERE 0j « (I - 1)w/2N 

(c) INTC=2 

1 and 2 in the A.P.R. 


32 





Figure 16 illustrates the use of NODE CARDS to fora 
chordwise regions on a section for a wing patch. Figure 16 (a) is 
for the case of a simple faction with two chordwise regions while 
Figure 16(b) shows a m o re c pm plicated case with five chordwise 
regions. This case includes an illustration of the more general 
copying capability to be described in 4. 5. 7.1. 

The chordwise region option can also be used on sections 
generated internally by the code under INMODE values of 5 through 
8. For example, INMODE = 5 causes a NACA four-digit equation to 
generate a symmetrical section. (Note: the four-digit equation 

has been modified in the quadratic coefficient to give zero 
thickness at the trailing edge: 

z i = + TC* ( 1.4845 yJTi - .63xi - 1 . 758xi 2 
+ 1 .4215xi 3 - . 518xi 4 > • 

The basic points on the section are generated with 

xi = 1.0- s i n 8 j ; 0 < 9i — 71 * 

(Note: this determines panel spacing if NPC = 0 on CARD 14.) 

e = u - nve-'NiNT + e a 

a 

where ^B is the interval in 0 over the chordwise region and NINT 

is the number of intervals to be generated in that chordwise 

region. (Warning: if NPC = 0 on CARD 14, then MINT is the 

number of panels; default - 70.) 

is the value of 0 at the start of the chordwise region. 
The 0 interval, Ve, over the region is evaluated as: 

V0 = sin~l ( 1 - | XRB I ) ~ 0 

a 

The user supplies the values for HINT and XRB on CARD 13. XRB is 
the x-station for the end of the required chordwise region; it 
takes a negative sign if the region ends on the lower surface. 

Figure 17 illustrates the card set required for the INMODE = 

5 option using four chordwise regions. It also covers the case 

for INMODE = 7 which generates a biconvex section. 

Figure 18 illustrates the case for IN MODE = 6. This 

generates a semi-ellipse with a unit horizontal semi-axis in the 
S.C.S. A similar chordwise region option is available as in the 
IN MODE = 5 case and uses the CARD 13^14 information, see Figure 
17 . 
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IF A SET OF POINTS IS DEFINED IN 
ANOTHER COORDINATE SYSTEM, USE THE 
"MOVE" PARAMETER ON CARD 14 AND 
SPECIFY TRANSFORMATION ON CARD 14B 


original PA ! 

m POOR QU, 


HOLD THESE AS PANEL CORNER POINTS 
(NPC“0) TO MATCH LOWER-SURFACE 
POINTS IN THIS REGION 



IF POINTS ALREADY INPUT FOR AN EARLIER SECTION 
USE A NEGATIVE NODE VALUE ON CARD 14 AND 
IDENTIFY THE PATCH, THE SECTION AND THE STRING 
OF POINTS IB | LB ON CARD 14A 



CARD llf-ft 

(BECAUSE NODEC. NEGATIVE 
ON CARD 14) 


Figure 16. Concluded. 

(b) Complicated Section. 
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HORIZONTAL SEMI AXIS = 1.0 

VERTICAL SEMIAXIS = H (INPUT ON CARD 12(F)) 


SIN 0j 


z, = H cos dj 


0, - d A + (i - l)Ad 

Ad * (0 B - 0 a )/nint ) 

dg = SIN -1 j XRB | 


NINT AND XRB 
INPUT ON CARD 13 


dg = V - dg IF XRB NEGATIVE 


d A , d fi ARE THE VALUES OF d AT THE 

i BEGINNING AND END/ RESPECTIVELY/ 
OF A "CHORDWISE" REGION. 


Figure 18. Illustration of Semi-Ellipse Option, INMODE=6. 
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A complete ellipse can be generated as a section using the 
I HI; ODE = C option. Figure 19. In this case the complete 
horizontal axis stretches from 0 to 1 . 0 in the S.C.S. The 
vertical semi - axis, II, is supplied by the user on C AS. D 12(f). 
The chordwise region option using CARDS 1 3 -"1 4 is available in a 
similar form to that shown in Figure 16. 


4.5.6 Spanwis e Regions 

Sections defined within each patch may be assembled in a 
number of spanwise regions for the purpose of controlling panel 
density and spacing in the spanwise direction. In forming 
spanwise regions, sections defined by the user take on a similar 
role to that of basic points in the chordwise regions. Although 
the options available for the spanwise regions are essentially 
the same as described for the chordwise regions in 4. 5. 5. 3, the 
two are applied completely independently; for example, the user 
may request automatic paneling in the chordwise direction and 
manual in the spanwise direction. As in the case of chordwise 
regions, spanwise regions are used only as an input convenience 
and are discarded once the paneling is complete. 

Spanwise regions must end at user-defined sections, called 
N ODE SECTIONS, Figure 20. These usually coincide with kinks in 
the spanwise direction on the patch planform, but the user can 
place one whenever he wishes to change the form of the paneling 
or to change between manual and automatic paneling in the 
spanwise direction. For convenience, the spanwise node 
information is included on the section card (CARD 11) together 
with the section transformation information (4. 5. 5.1). The 
function of the spanwise region node quantities, NODES, UPS, 
I NTS, follows closely the description in 4. 5. 5. 3. MODES, 
however, must be set to zero (b\Jank) on the first section of a 
patch and on all intermediate input sections that are not node 
sections. (NPS and I NTS are then inactive.) The last section on 
a patch is identified by a NODES value of 3, 4 or 5; 4 is used if 
the patch is the last one on a component and 5 is used if the 
patch is the one on the configuration, in which case the present 
section completes the basic description of surface geometry. 

An exception to the above occurs in the case of a body of 
revolution. To execute this option a negative sign is placed on 
the MODES value on the first section card. The section basic 
points are then generated in the S.C.S. z,x plane. Figure 21. 
the section is rotated about the x-axis through an angle given on 
CARD 15; this generates the 'spanwise' paneling in accordance 
with the NPS, INTS data on the section card. Multiple spanwise 
regions may be used, if required, as illustrated in Figure 21. 

The total number of panels defined (manually or 
automatically) across each patch in the spanwise direction is 
monitored by the program ana becomes the number of panel columns, 
NCOL, for that patch. In view of the ease of generating panels. 
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Xj * 0.5 (1.0 + cos 0j) 

Zj = -H SIN 0J 

0j “ + d - 1)A0 

A 0 = (0 R - 0J/NINT ) 

( XRB, NINT INPUT , 

0 B * cos' 1 ( <J XRB | - .5)/. 5) ) 0N CARD 13 

fig - 2TT- 0g IF XRB NEGATIVE 


Figure 19. Illustration of Full Ellipse Option (INMODE=8) . 



INTERMEDIATE INPUT SECTION 



■0 



THE SECTION DEFINITION 


Figure 21. Body of Revolution Option. 
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the code also monitors the running total of panels, and if a 
limit is exceeded, the calculation terminates with an appropriate 
error message* The limit is set internally by the storage 
capacity, but the user is given the opportunity (with variable 
NPNMAX on CARD 3) to override that value with his own estimate of 
the total he iajtends. to use for that case. In the event of an 
input error, this will avoid the inadvertent and expensive use 
of, say, 1,000 panels when the user intended using only 100. 


4.5.7 S^e i c 1 a 1 ^Routines 

The geometry routines described abov^e may be applied for all 
patches on the configuration surface; however, special routines 
have been provided to reduce user input and, in particular, to 
avoid duplicating information already supplied. These routines, 
which are described below, are optional. 

4 . 5 . 7 . 1 Copying Rout ine 

We have already seen (4. 5. 5. 2) a copying facility accessible 
at the section input level. This copies over a complete section, 
including the chordwise region information, and has, therefore, a 
rather limited application. More general copying routines are 
provided and are activated at the basic point level to copy 
STRINGS OF BASIC POINTS, rather than complete sections. This 
capability allows a new section to be assembled fro m £ar t s o f 
previously defined sections. Several strings of basic points may 
be assembled from a number of previously defined sections and the 
points selected need not follow the same direction as originally 
specified. Furthermore, the copied strings of points may be 
intermixed with strings of manually input basic points to 
complete the new section. 

For this copying inode, the value of IN MODE on the section 
card 11 must be in the range 1 to 4. The copying is activated by 
inserting a MODE CARD having a NEGATIVE sign on NO DEC. This is 
regarded as a DUIiliY node card because it does not necessarily 
terminate a chordwise region (see below). The negative value for 
NGDEC determines the action at the end of the copied string of 
basic points. If NGDEC = -1 or -2, then the last copied point 
becomes the end of a chordwise region on the new section and 
signifies, respectively, continuous or discontinuous slope onto 
the next chordwise region. We then continue to specify further 
basic points, or, by inputting another negative node card, we can 
copy another string of basic points, and so on. If N0DEC = -3, 
then the last copied point in the string completes the new 
section. 

If the user does not require a chordwise region to end at 
the last point in a copied string, then he sets 1 1 0 D E C = - 4 when 
he initiates the copy. V/ hen the string has been copied over, the 
program then expects to receive further basic points to complete 
the chordwise region or another negative node card can be used to 
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copy inothe r string of points, and so on. Clearly, if NODEC = 
”4, then the NPC and INTC values on that NODE CARD are inactive 
and aay be left blank. 

Whenever a negative node card is inserted, it most always be 

followed by a COPY CARD (CARD 14A) conta ining the following 
inf oraat ion (four integers) de f ining the location of the required 
string of points, IPCH, I SEC, IB, LB . 

IPCB is the patch number containing the required points . 

ISEC is the section number relative to the start of that 
patch. 

IB, LB are, respectively, the first and last basic point 
numbers ( inc lus ive ) defining the string . The 

numbering is relative to the start of the section 

ISEC. 

(Note : these subscripts can be obtained during an MSTOP = 1 run 

(CARD 2 ) with IPROGOM print control set on CARD 2A. 

Thus, even in a compl icated conf igurat ion, it is relatively 
easy to specify a s tr ing of basic points. 

This option offers not only an alternative to the earlier 
copying routine, but also a more general capability because the 
copying is initiated a t the basic point input level, rather than 
a t the section input level. For example, the comple te copied 
section need form only a part of the new section, it being 
possible t o have other basic points, both before and after the 
copied string. In addition to this, the ability to break the 
copying into strings of points allows a new distribution of 
chordw i se regions and pane 1 ing to be selected. 

One restriction must be considered when using this copy 
routine — the new section's value for INPUT muz t coincide with the 
INPUT values on sections from which strings of points are to be 
copied. This restriction ha s not posed a problem so far, but if 
it doe s , it would not be too difficult to remove it. 

This copy routine is used in the w ing -body sample case in 
the Appendix. The fuselage lines a t the w ing-body junction are 
copied from the first section on the wing patch. 

4 . 5 . 7 . 2 Automatic Patch Generator 

Patches covering tip edges, flap edge, cutout s, etc. can be 
input by the user a s ordinary patches, but this can get tedious. 
Optional automatic procedures have ben instal led which simplify 
this input by generating a complete patch within the code . This 
AUTOMATIC PATCH GENERATOR, or A.P.G., is initiated at the patch 
input level by inserting a non-zero value for parameter, MAKE, on 
the patch CARD 10. The value of MAKE ident if ie s the patch number 
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on the edge of which « closing patch is to he generated. The 
sign of MAKE determines whether t he new patch is on side 3 
(positive) or side 1 (negative) of the basic patch. The patch 
referred to should be a folded patch although sides 2 and 4 need 
not necessarily meet for the A.P.G. More importantly, an equal 
number of points (or panels if NPS = 0) should have been used on 
the upper and lower sides of the edge section. 

Cons ider, for example, a tip-edge patch. Here we have 
already defined the patch repr e sent ing the main surface. The 
end section of that patch provide s the BASE SECTION from which 
the A.P.G. creates the new patch. Figure 22(a), according to user 
instructions. When the A.P.G. has been activated, the next card 
( CARD 16) must contain the following: 

NPC , INTC, KURV, NPTIP, NPS, INTS 

Referring to 4.5.5 and 4.5.6, the generated patch has one 

chor dw i se region with NPC panels spaced according to the va lue of 

INTC. 1 1 has one spanw i se region with NPS panels spaced 
according to INTS. In this case, the basic points ( in the 

C.C.S.) defining the base section are brought over. I f NPS = 0 , 

then the main program uses the actual pane 1 corner points on the 
base section; clearly, this gives an exact ma tching of panels 
across the patch edge. The values of NODES must be either 3 , 4 
or 5, depending on the location of the patch in the input. The 
function of KURV is de scribed below. 

Sections def ining the new patch are created automatically 
from the base section coordinates. The contour of each section 
generated may be either a straight 1 ine ('square-cut' tip) or an 
ellipse or triangle, depending on the value of the quant i ty, 
KURV, supplied by the user. I f KURV is 0 , sets of basic points 
are generated on straight lines joining upper and lower points on 
the base section. The same number of points is created even if 
the interval across the base section is zero ( e . g . , a t the 
leading and trailing edges). Figure 22(b). 

I f KURV is 1 , the basic points are created on semicirles 
having diameter e qua 1 to the local 'thickness' of the base 
section. 

If KURV is 2, the basic points are created on semi-ellipses; 
the base section local thi ckne s s provides one axis, while the 
horizontal semi-axis is derived from additional user input , Fig- 
ure 22(c). A pi an f or m shape is input using a set of coordinate s, 
s * , y * , i = 1, NPTIP, def ine d in a convenient local coordinate 
system with origin at the trailing edge. Figure 22(d). The scale 
and point distribution are compl e te ly arb itrary, so the points 
may be conveniently measured from a pi anf orm view of the wing. 
The program scales the shape to fit the length of the basic 
section and interpolates to find the local sem i-axi s for each 
ellipse. 
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IF. MAKE * -NPCH ON CARD 10 



IF MAKE = +NPCH ON CARD 10 
TIP PATCH GENERATED 

(A) PATCH LOCATION FROM SIDE 3 


STRAIGHT-LINE SECTIONS GENERATED 
BETWEEN UPPER AND LOWER POINTS. 
PANELING IN THIS DIRECTION IS 
DETERMINED BY NPC, INTC ON 
CARD 16 ( 1 . E . , "CHORDWISE") 


y GENERATED PATCH 

* 

SIDE" i 
(POINT) 

(B) FLAT PATCH OPTION (KURV=0) 


Figure 22. Automatic Patch Generator (A.P.G.). 

(a) Patch Location. 

(b) Flat Patch Option (KURV=0) 
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KURV * Is SEMI-CIRCULAR SECTIONS 
GENERATED BETWEEN 
UPPER AND LOWER POINTS 


KURV = 2: SEMI -ELLIPTICAL SECTIONS 
GENERATED. USER MUST 
SUPPLY PLANFORM SHAPE 
(s r Yjj I * 1, NPTIP 

ON CARD 16A) 
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Figure 22. Concluded. 

(c) Semi-Circular and Semi-Ellipse Sections 
(KURV=1 and 2) 

(d) Triangular Sections (KURV=3) 
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A KURV value of 3 is treated in a similar way to the seiu- 
ellipse case, except the tip planform description must now 
include vertical offsets which locate the apex of each triangular 
section (see Figure 22(d))® 


4 . 6 P a.n e 1 

4.6.1 Pane 1 Gene r a t j. on 

When the basic geometry has been specified, the panel corner 
points are assembled, patch by patch. A temporary set of 
chotdwi.se points corresponding to panel corners is first 
assembled on each of the defined sections. This is performed in 
each chordwise region in turn ( 4.5 .5 .3 ) and interpolation is used 
when the A.P.R. is requested ( i . e . , when NPC > 0). The form of 
the interpolation used by the A.P.R. depends on the number of 
basic points available in the chordwise region, including the two 
end points. The code augments this number by talcing a basic 
point from a neighboring chordwise region if continuous slope has 
been specified onto that region (i.e., NO DEC = 1 on this or the 
previous region). The A.P.R. takes the available set of basic 
points and first eliminates zero length intervals, then, 
depending on the number of basic points left, i.e., one, two, 
three or more, it uses, respectively, constant, linear, quadratic 
or biquadratic interpolation to generate the panel points. 

Vhen the temporary set of chordwise points is complete for 
all sections on a patch, corresponding points on each section are 
joined by lines called SPANt/ 1 SE GENERATORS. The panel points 
along each spanwise generator are then assembled in a similar way 
to that described for the chordwise direction, but now based on 
the spanwise region information. The interpolation routine is 
now applied along each spanwise generator in each spanwise region 
where the A.P.R. lias been selected. The new set of (spanwise) 
points are actual panel corner points from which the panel geome- 
try is generated. 

The fact that we input just one set of spanwise region 
information for a patch means that the same spanwise 
interpolation format is used on all the spanwise generators on 
that patch. Thus the A.P.R. in the spanwise direction has lost 
generality compared with the chordwise capability; however, this 
loss is not serious (and to avoid it would require considerably 
more input). 
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4.6.2 Panel Geomet ry 

The four corner points, Rj , i = 1 , 4, specifying a panel 
quadrilateral are assembled in the same sequence as the corners 
on the parent patch. Figure 23. The panel's control point, R c 
(where the boundary condition is applied) is the near, of these 
four points. 

The parallelogram formed by joining the mid- points of the 
sides in^sequence provides the near, plane of the g^anel. The unit 
vector, n, normal to this plane and unit vectors, 1 and at, within 
the plane form the panel's right-handed orthogonal unit vector 
system for 1 o c_g 1 coordinates. This system takes the panel 
control point, P. c, as origin and the vector, % , passes through 
the midpoint of side 3. 

The four corner points are projected onto the mean plane to 
form the flat panel used for the velocity potential influence 
coefficients. The distance between the actual corner points and 
the mean plane (this is the same for all four corners) is a 
measure of the amount the panel is skewed. A warning is printed 
when this value becomes large. 

The program forms the half-distance. Slip, S II Q from the 
control point to ; t h e midpoints of sides 2 and 3, respectively (or 
sides 4 and 1), see Figure 23. These distances are used both as 
a measure of panel size (SUP + SNG) and also as a measure of 
surface distance in the evaluation of doublet gradients. Within 
a patch the side midpoints of neighboring panels coincide and so 
linked S U p values and linked S 11 Q values give a close 
approximation to surface distances between panel control points. 


4.6.3 Pane l N eighbors 

The program keeps two files on panel neighbor information. 
This information includes the subscripts N A D 0 R * ; i = l,4 (in the 
sane sequence as the sides) for each panel's set of immediate 
neighbors and NABSI 1) i , the sides of those neighbors adjacent to 
the panel, see Figure 24. Initially, for panels within a patch 
this information is rather trivial; only the information across 
patch edges is significant. The latter information is formed by 
an automatic procedure performing a neighbor search. All this 
initial information is stored in one file while a second file is 
formed in which the neighbor information is modified at expected 
jumps in the doublet distribution. The files are referred to as 
'before wake shedding' and 'after wake shedding', respectively. 
Both sets can be printed using the IPRNAP. option on CARD 2A. 
Temporary cuts in the neighbor relationships are indicated in the 
second file by a zero NABOB value and a negative NAB SID. A panel 
side on a plane of symmetry has a negative NABSID and its own 
subscript in the NABOB value in both neighbor files. 
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Figure 23. Panel Geometry. 






The panel neighbor information can be codified by using the 
options on CARD 8. The new information override s the automatic 
neighbor information. This facility is used to correct S.LL2.LS. 
made by the automatic procedure in awkward patch junctions or to 
obtain special effects. I t should be noted that the panel 
neighbor information does not normally affect the doublet 
solution--it mainly affects the doublet gradients evaluated in 
the analysis procedure. An exception to this would occur if 
neighbor errors occurred along a wake-shedding line--this could 
then cause a doublet solution error because wake shedding panel 
information (see later in Figure 27) would be wrong. 


4 . 7 Wake 

4.7.1 Wake-Grid -Planes 

The wake is formed after all the surface patches have been 
panelled and the panel neighbor information stored (before wake 
shedding). Y/ithin the code, the wake geometry is described in a 
set of cross-flow planes called WAKE-GRID-PLANES, whose x- 
stations are prescribed by the user. Figure 25. Spacing between 
w ake-gr id-pl ane s should be small where intense vortex roll-up and 
vor t e x^*s ur f a c e interaction is anticipated. The intervals between 
wake-grid-planes can be progressively widened with distance 
downstream beyond the region of interest. 


4.7.2 Wake Parameter s 

At the input stage the wake has the following parameters: 

Vi'N AliE , IDENTW , IFLEXW, IDEFW 

WNAME is a text identification for user convenience, while IDENTW 
is the wake type. Two 'wake types' are considered in the present 
document: IDENTYV = 1 is a regular wing-type wake with constant 

doublet distribution down each strip in the streamwise direction; 
and IDENTW =4 is a jet-type wake whose stripwise doublet 
distributions are linear in the streamwise direction. The 
doublet gradient is specified as the jump in tangential velocity 
across the type 4 wake sheet. the IFLEXY7 parameter is used to 
specify relaxable (0) or rigid (1) wakes. The I DEFY? parameter 
determines the way the wake FIXED BASELINE is to be defined. 
IDEFW® 0 requires strings of WAKE-SHEDDING P ANEL S-- d e s c r i b e d 
below — while IBEFW = 1 requires a set of separation points. The 
latter option is being checked out at this time. 
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Inside tie program wakes have additional parameters 
associated primarily with WAKE COLUMNS of WAKE PANELS. Figure 26 
shows the arrangement of these in relation to the wake-grid- 
planes. The view is from fi.l> 5 >v£; i . e . , wake panel normals are 
directed towards us. The first wake panel ( a t the top) on each 
wake column (NWC) has subscript IWPAN(NWC). The total number of 
columns on a wake is NWCOL and on the whole configuration. 
NWCOLT. 

Each wake column is associated with a set of four surface 
panels, two upstream (or Upper) and two downstream (or Lower) of 
the separation line. These are KWPU, KV/PUU, KWPL, KUPLL, Figure 
27. The corresponding surface distances between control points 
and separation line are also stored SU, DSU, SL, DSL. The onset 
flow potential at these panels are also stored PI1IU, PHIUU. PHIL, 
PI1ILL . 


Each WAKE LINE (along an edge of a wake column) keeps the 
first wake-grid-plane, IWGP, intersected downstream of the 
separation point. The array, L5EQU, keeps the subscript sequence 
of wake lines intersecting the. wake-grid-planes. A negative sign 
is placed on subscripts of wake lines at the edges of each wake. 


4.7.3 Wake-Sh e dding Panels 

To locate the initial wake, the user identifies strings of 
wake-shedding panels, the side geometry of which defines the 
fixed base line of each wake. The basic information required to 
define a string of wake-shedding panels is 

ICWPACII, KV/SIDE , KV/LINE, KWPAN1 , KWPAN2 

KWPACH is the patch number and KWSIBE is the side of that 
patch that is parallel to and has the £aS£_Ml£®ii£S as the 
separation line on that patch. Figure 28. K W' LINE determines 
which line of panels (row or column) is providing the separation 
line within the patch. KW’LINE may be left blank ( i . e . , 0) if 
separation is from the patch edge, KWSIDE; otherwise it refers to 
the row number of column number (as defined in Figure 28) of the 
panels shedding the wake. The wake-shedding panels are regarded 
as being 'upstream* of the separation line; they are always on 
the ■Le.f.t, when looking in the direction of the separation line 
(Figure 28). 
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Figure 26. Wake Arrangement, 
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Figure 27 . Wake-Shedding Parameters . 
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Figure 28. Separation Line Definition, 
(a) Wing. 









SIDE li 


COLUMNS 


SS* ® 


SIDE 1 


© 


ROWS (g) 


© ® ® ® 


a 

B 

B 

B 

B 

B 

m 

B 

B 

HI 

m 

MS 

B 

B 

□ 

m 

m 

m 

a* 

a 




rt 

st 

ti 

m 



m 

m 

m 

St 

Bi 

it 

m 

a 

a 

B 

V* 

B 

7 m 

______ _ 


m 


PATCH "KWPACH" 


SIDE 2 


SIDE 3 


SIDE "DIRECTIONS" 


SEPARATION LINE "DIRECTION" 
, WITH SEPARATION PANELS ON 
"LEFT" SIDE 


The separation line shown above is described by 7 sets 
of wake-shedding panels (shaded) as follows: 


Set No. 

KWPACH 

KWSIDE 

KWLINE 

KWPAN1 

KWPAN2 

(Remarks) 

1 

1 

2 

2 ( row ) 

1 

2 

(Counting along 
row 2 in direction 
of side 2 

2 

1 

1 

3 (col) 

3 

4 

(Counting down Col 
3 in direction of 
side 1) 

3 

1 

2 

4 (row) 

3 

5 


4 

1 

1 

6 (col) 

5 

5 


5 

1 

2 

5 (row) 

6 

9 


6 

1 

3 

9 (col) 

1 

4 

(Counting along 
side 3 

7 

1 

4 

2 (row) 

1 

3 





Local counting, not absolute 
panel subscripts 


Figure 28. Concluded. 

(b) General Case. 
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KWFAN1, KWPAH2 are the first and last panels, respectively, 
in the present string of wake shedding panels. These are l.oc.a.,1 
numbers starting with 1 at the patch edge and proceeding along 
the row or column i. u_t. h£_d,i r e £.t .i o n_o.f _Jt h£„,l i.n£_b£i. ng^d e .£ i. n£d : 
using local numbers along the line avoids the complicated task of 
identifying a string of a.b s.o l,u jte. panel numbers which are not 
necessarily on a continuous sequence. For example, on patch 2 in 
Figure 28(a), the separation line is along row 3 ( K V LIME) 
parallel to side 4 and so the local panel numbering sequence 
starts with 1 at side 3 and proceeds to 5 at side 1; in this 
local sequence, the separation panels, KWPAN1, KWPAN2 are 
therefore 4 and 5, respectively. (The corresponding absolute 
panel subscripts are worked out inside the code.) A more general 
separation line path is shown in Figure 28(b); this illustrates 
most of the possibilities of this scheme for describing the 
separation line location. 

A separate string of wake- shedding panels must be defined 
for each patch crossed by the separation line. 


4.7.4 S t reamw i s e Wake Line s 

At the beginning of each string of wake- shedding panels and 
at the end of the last string on each wake, the user must specify 
the geometry of a STREAMWISE WAKE LINE using STREAM? I SE WAKE 
POINTS (SWPX, SVPY, SWPZ), Figure 29. The function of these 
points is similar to that of basic points defining patch sections 
(4. 5. 5.2). MODE cards are used here also and allow the user to 
augment his input points with interpolated points. As a minimum, 
one streamwise wake point must be provided at some station beyond 
the last wake-grid-plane. (The upstream end is taken 
automatically from the fixed baseline identified by the string of 
wake-shedding panels.) 

The streamwise wake points are defined in a local coordinate 
system parallel to the global coordinate system but with origin 
at the local separation point on the fixed baseline. The program 
relocates the streamwise wake lines into the global coordinate 
system. 


An option is provided to copy and rotate the previous 
streamwise wake line. The rotation (DTH) is about the x-axis of 
the local system and is in the same sense as the body of 
revolution option. Figure 21. 



FIRST POINT ON EACH WAKE LINE 
IS TAKEN AUTOMATICALLY FROM 
THE FIXED BASELINE 
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4.7.5 Wake Panel Ppin_i_s 

When the basic wake information has been supplied, the 
program generates streamw i se sets of coordinate s describing the 
path of each of the defined wake lines according to the user's 
instructions on the node cards. Biquadratic interpolation is 
used to generate these points based on the user supplied 
coordinates. Linear interpolation is used in the span wise 
direction to generate intermediate wake lines between the defined 
lines. Wake lines are formed corresponding to each of the wake- 
shedding panels. As each wake line is formed, its intersections 
with the WAKE GRID PLANES are computed and these intersection 
points become the wake panel corner points. LINEAR interpolation 
is used in the computation of these intersection points. The 
user should bear this in mind when selecting the number and 
location of both basic points and streamwise lines. It must be 
emphasized, however, that this information is used only to define 
the preliminary wake for the purpose of the first solution-- 
thereaftcr, the wake relaxation routine will redefine the wake 
geometry at each iteration, In the unsteady mode, this basic 
wake is transported downstream as new panels are formed at the 
separation line. 

The points generated in the w ake-gr id-pl ane s are the actual 
wake panel corner points. Figure 26. These are treated in 
exactly the same was as surface panels (see Section 4.6). 



5,0 OFF-DODY VELOCITY SCANS 


5 , 1 O ptions 

Of f -body velocity calculations are performed at user 
selected points. Normally, the points are assembled along 
straight scan lines. These lines in turn may be assembled in 
planes and the planes in volumes. The shape of the scan volume 
is determined by the parameter MOLD. Two options are currently 
available : 

MOLL-1 Allows a single point, points along a straight 

line, straight lines within a parallelogram or 
within a parallelepiped. 

MOLD =2 Allows points along radial lines in a cylindri- 

cal volume. 

A number of scan boxes can be specified, one after the 
other. A M0LD=0 value terminates the off-body velocity scan. 

The location of scan lines within the scan boxes can be 
controlled in the input. The default option is equal spacing 
generated along the sides of the box. The location of points 
along a scan line can be controlled also in the input, but again, 
the default option is equal spacing. 

If a scan line intersects the surface of the configuration, 
an intersection routine locates the points of entry and exit 
through the surface paneling. The nearest scan line points to 
each intersection are then moved to coincide with the surface. 
Points falling within the configuration volume are identified by 
the routine to avoid unnecessary velocity computations. These 
points are flagged in the printout. The input parameter MEET 
controls the calls to the intersection routine: I!EET=0 (default) 

makes the routine active while IIEET-1 switches it off. Clearly, 
if the velocity scan volume selected is known to be outside the 
configuration volume then switching the intersection routine off 
avoids a lot of wasted computation. 

The velocity routine VEL computes the velocity vector at 
each point by summing the contributions from all surface panels, 
wake panels, image panels (if present), onset flow, etc. The 
routine includes near-field procedures for dealing with points 
that are close to the surface panels or wake panels. The near- 
field routine involves a lot of extra computation, if, therefore, 
it is known that a scan box is well clear of the configuration 
surface and wake then it is worthwhile turning off the near-field 
routine. This is controlled by the input parameter NEAR; NEAR=0 
(default) keeps the near-field routine active while HE AR= 1 
switches it off. 
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Clearly, because of the MEET and NEAR options, it may be an 
advantage to divide large scan volumes into a number of smaller 
bones, many of which could have the near-field and intersection 
routines switched off. 

The velocity scan data (X, Y, Z, VX, VY. VZ, V, CD for each 
point on a scan line is written away to a plot file (TAPE 7) 
together with other data. Since this provides a ready means of 
examining the off-body data quickly and since a large amount of 
printout can be generated by the off-body scan routine, then an 
option is provided to just print out a sample of the calculated 
data. The frequency of this print by plane, line or point is 
specified by the input parameters, INCPF.I, INCP11J, IIICPI.K, 
respectively, A zero (default) causes all the appropriate data 
to be printed while a value of N , say, causes the results to be 
printed every Nth plane, line or point as appropriate. 


5 . 2 M0LD =1 . Sk ewed Box 

The M 0 L B = 1 option allows a very general scan volume to be 
defined, Basically, three edge vectors are input and a 
parallelepiped is constructed. Figure 30(a). However, the 
position vectors defining the edge vectors are input one at a 
time giving the option to stop the input after each point, thus 
allowing a single point, a straight line or a plane to be 
considered in addition to the volume option. 

Referring to Figure 30(a), the coordinates of the first 
corner point are input 


xo, yo, ::e, i:p 

NP controls the various options as follows: 

MP=1 This is a onc-p o int box and therefore has no 

further input (but it can be followed by another 
scan box with any HOLD options) 

KP>1 A line (^1) is generated after inputting the 

second corner point: 

XI, Yl, 21, NP 1 

where HP1 is the number of points (equally 
spaced) along the side SI 

NP=2 This is a line scan only and so there is no 

further input for this box 
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IN THIS BOX ARE PARALLEL 

to this Line) 


Figure 30. Velocity Scan Volumes. 


(a) Skewed Box. 



i:P > 2 The third corner point is input : 

X2, Y2, Z2, NP2 

4 

TIi i s defines side S2 witli NP2 points along it. 

HP=3 This is a scan plane only and so there is no fur- 

ther input for this box. The scan plane is _£he 
parallelogram generated from the two sides, SI, 

S2. Eacli. scan line within the plane is parallel 
to side SI* has the same point distribution as SI 
and is located by one of the points on "&2. 

!IP = 4 The fourth and last corner point is input: 

X3 , Y 3 , Z3 , 1JP3 

This defines side si* with NP3 points along it. 

These points locate a set of parallelograms with 
sides parallel to "Si, "§2 and each one having a set 
of straight scan lines as described for IiP-3. 

The points along each scan line are, by default, equally 
spaced. However, negative signs can be placed on 'K.P1, HP2 or NP3 
and then the point locations along the corresponding sidcs(s) 

must be input; e . g . , C a i , i = l, I Ill’ll, where 0 < a* < 1.0, locates 
points along vector SI. 

5 . 3 H0LD = 2 j Cy lj,n d r j. c.a.1 Voluuc 

The cylindrical volume scan requires a minimum of three 
input cards. The first card defines the location of the center 
of a circle (XI, Y 1 , 21); the inner and outer radius there (HI, 

RG); and the first and last 0 values (01# 02^' measured from the 

vertical. Figure 30(b). 

The second card defines a second circle plane with center at 

X2, Y2, 22 , respectively, and with inner and outer radii, Rl2' 
# respectively. The values of 6i , ©2 are assumed the same as 

on the first card, 

The last card specifies the number of points (UAL) along the 
cylinder axis between X 1 , Y 1 , Z 1 to X2, Y 2 , 22; the number of 0 
values (NTilETA) between $2# and the number of points (NIL AD) 

along radial lines between local inner and outer radius. The 
points are located using equal spacing by default. However, as 
before, a negative sign placed on the number of points allows 
these points to be specified in an arbitrary distribution. 

A point on a scan line in the cylindrical volume is located 
as follows (referring to Figure 30(b)). 
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First define vector D1 along the cylinder axis, i . e . , 

D1 = tl2 - RI 


-4 

From this the unit vector, h, is formed, also the unit vector, n, 
which is in the horizontal plane and which is normal to the 
vertical plane containing h, i.e.. 


n = k_2_h 

It - ii 

^ A" 

Also, the unit vector, g = n A + h, is + generated forming a local 
unit orthogonal vector system, n, h, g. 


The unit radius vector, u, at angle 0 to the vector g is 
given by: 


or 


A t s t sin 6 

4 *f _ 4 , 

u = g co s 0 + n s 1 n 0 


Thus, the position vector of the Kth point on the J'th radius 
vector at the I** 1 location along the cylinder axis becotnes 

t » ^1 + a R (K) * du(J) K= 1 , NEAR . 

where 

tl =tl + a( I ) * t>l + III(I) * u(J) J=l, NTIIETA 

1=1. NAL 


u(J) is evaluated with 


e(j) - e 1 + a T( j) * (6l _ e2 ) 

HI (I) = RI i + a( I ) * (Rig - RIi ) 


and d = EOj 

is the length of the 


- Eli + ct(I) * (R0 2 - RI 2 ~ ROi + RI i ) 
present radial scan line. 
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6.0 


RUNNING THE PROGRAM 


01! POOR QUALITY 


$ . 1 Job Control Card* 

A typical J CL file for running the VSAERO program on the 
CRAY computer at NASA Ames Research Center is listed below : 


JOHNDOE, MICE. 

USER, XXXXX , XXX . 

JOB. JN* JOHNDOE , T»200, US-XXXX. 
ACCOUNT, AC-XXXXX. 

ACCESS, DN«$BLD, PDN-VSLGO, ID=XXXX. 
ASSIGN, DN«=DATA, A-FT05. 

LDR. 

REWIND, DN-FT06. 

SAVE, DN=FT06 , PDN-WINGDO, ID-XXXX. 


The binary program is stored here as a direct access file 
called VSLGO. The input data file is called WING and the output 
file is saved as a direct access file, W INGDO. 
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6 . 2 Plotting Data 


OF 


la order to plot the date interactively the plot file (TAPE 
7 ) amt be saved after the ran. To do this requires the 
following sdditionel cards in the J CL: 

REWIND, DN-FT07. 

SAVE, DN-FT07, PDN-WINDDP, ID-XXXXX. 

This set of cards should be located after the LDR card in 
the basic JCL, see 6.1. The plot file is then saved as a direct 
access file called 1IN6DP, which can be plotted using the routine 
OMNIPLT. 

The OMNIPLT rout ine steers the user through plots of panel 
and wake geoae try, sections through the surface pressure and 
velocity distributions, on-body stress 1 ine s (if present) and 
boundary layer data (if present), see 8.2. 

A typical day file follows . 
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CRAY-X SERIAL, 34 "ADVANCED COMPUTATIONAL FACILITY 06/17/82 

'V fP? .. • ; ' ’ , 

CRAY-1 OPERATING SYSTEM COS 1 10 ASSEMBLY DATE 06/09/82 


JOB, JN=BR I AN. T-200/ US* XXXXX 
ACCOUNT. 

ACCESS, DN«*2LD. PI>N*VSL.GOMAY, ID* XXXXX 

PD001 - ACCESS VSlGQMAY ED«0001 COMPLETE 

LOR. 

LD010 - BEGIN EXECUTION 

FTOoJ - STOP IN MVP 

REWIND. DN*FT07. 

SAVE, DN*FT07, PDN*WN©PLT, ID* XXXXX 

PD001 - SAVE WNGPLT ED*0009 COMPLETE 

END* OF JOB 


UOBNAME BRIAN USER NUMBER XXXXX 

TIME EXECUTING IN CPU - 00:02:38.2016 
TIME WAITING TO EXECUTE 00. 00.43. 3733 
TIME; WAITING FOR I/O — 00. 01: 01. 3*141 
MEMORY USAGE - — 57. 73944 MWDS-SEC 

: BLOCKS MOVED 32083 

PHYSICAL I/O REQUESTS - 7366 

TOTAL COST INCURREO — « . 32.29 




ORIGINAL PAGE IS 
OF POOR QUALITY 

7 . 0 INPUT DESCRIPTION 
7 . 1 Innut Summary 

The input is divided into the following parts: 

(i) EASIC INPUT 

General information, operating mode, onset 
flow, reference conditions, special options 

<ii) PATCH GEOMETRY 

Description of configuration surface in com- 
ponents, patches, sections, basic points, 
etc., for panel generation 

(iii) Y/AKE INPUT 

¥/ake-gr id-plane s , type of wake, wake separa- 
tion line, initial streamwise geometry 

(iv) SURFACE STREAMLINE INPUT 

Location of starting point for each surface 
streaml ine 

(v) BOUNDARY LAYER INPUT 

Reynold's number, etc. 

(iv) OFF-HODY STREAMLINE INPUT 

Location of starting point and required up- 
streamed owns tr earn distances for each off- 
body streaml ine 


In the following description, the input variables are first 
listed in 7.2 for each of the above parts. Then, 7.3 gives a 
detailed description of the function of each input variable. 
This is followed in 7.4 by an input flow chart to help with the 
assembly of the input data file. 
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1 • 1 iRRfi tJVar iable Lis t 


ORIGINAL PAGE IS 
OF POOR QUALITY 


SfeJLi.jg. -JLa |) ut Sammarv 


Card_No . Variables 

1 Text 

2 IPRI, IPRLEV, IPRESS, MSTOP, MSTART, MODIFY 

2A IPRGOtl, IPRNAB, IPR1VAK, IPRCPV, IPRPPI (only 
if IPRLEV=5 on CARD 2) 

3 MODE, NPNMAX, NRBLiAX, ITGSMX, IMERGE, NSUB, 

NSPMAX, NPCIIAX 

3 A NROV/B ( I ) , 1 = 1, llJRBKAXl (only if NRDMAX< 0 on 

CARD 3) 


Format. 

20A4 

615 

515 

315 

1615 


4(a) NY7IT , NVPI, ICLTYP (if M0DE=1 on CARD 3) 315 

or 

4(b) NT, NEC (if MODE =2 on CARD 3) 215 

4 A (only if NVPI>0 and IBLTYP=0 on CARD 4(a) 

(i) UPSETS 15 

(ii) NPCKEL , NBCOL, (KOMI), 1=1, NBCOL) 1615 

(Number of 4A(ii) cards = MPSETS) 

If MSTART >0 and MODIFY=0; this is the end of 
the basic data on a restart run. 


5 RSYK, RGPR, RMF , RFF, RCORE , SOLRES, TOL 7F10.0 

6 ALP.BG, YAWDEG, REACH , VMOD, COI1FAC 5F10.0 

6 A ALBAR, RFREQU, HX, MY, HZ (only if IiODE=2 on 5F10.0 

CARD 3) 

7 CBAR, SREF, SSPAN, RIIPX, RMPY , RMPZ 6F10.0 

8 DORSET , NVOHT, NPASBM, JETPAH, NBCI1GE 515 

8A ( HORPCII ( I ) , NORF(I), NORL(I), NOCF(I), 515, 

IIOCL(I), VNORK(I) , ADBB(I), 1 = 1. HORSET) 2F10.0 

(only if HORSET) 0 on CARD 8) 

8B ( i ) VOET i( o n 1 y if NV0RT>0 F10.0 

> on CARD 8) 

( i i ) ( RXV ( I ) , RYV(I), RZV(I),J 3F10.0 

1=1, HVORT+1) 
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Card No . 

Variables 


Foruat 

8C 

( NPSPC1I ( I ) , NPSRF(I) , NPSRL(I), 
NPSCL(I) , 1 = 1, NPASUII ) (only if 
CARD 8) 

NPSCF ( I ) , 
NPASUM > 0 on 

515 

8D 

( JETPCI! ( I ) , JETRF(I) , JETRL(I) , 
JETCL(I) , VIN(I), VOUT(I) . 1=1. 
(only if JETPAN>0 ON CARD 8) 

JETCF(I) . 
JETPAN) 

515 

2F10.0 

CO 

(EPAN(I), KSIDE(I), NEWHAB ( I ) , NEWSID(I), 1=1, 
NBCHGE) (only if NBCHGE >0 on CARD 8) 

415 



Patch Geometry Input, .Summary 


C a r d_N o . 

Variable 

Forma t 

9 

CTX , CTY , CTZ , SCAL. TIIET (component card) 

5F10.0 

9 A 

CPX, CPY , CPZ, CKX , CDY , CIIZ (only if SCAL<0 
on CARD 9) 

6F 1 0 . 0 

10 

I DENT, HARE. KOHP, KLASS, PNAIIE (patch card) 

415 , 
6A4 


(Note: If HAREM), go directly to CARD 16) 


11 

STX, STY, STZ , SCALE, ALF, TIIETA, INNODE, MODES 
NPS, INTS (section card) 

9 

12(a) 

DY, BZ. X ( INNGDE= 1 ) 

\ 


(b) 

EX, DZ, Y (IN1ICDE = 2)| 

3F10.0 

(c) 

BX, BY Z (INIiODE=3) 1 


(d) 

/ 

BX, BY, BZ (INMODE=4), 


( e ) 

TC, INPUT (IHHODE=5 or 7)1 

F10.0, 

(f) 

II, INPUT (INMODE=6 or 8) ) 

15 

(g> 

BX, RAD, TKET (INNODE=12) 

3 FI 0 . 0 

13 

XRE. HINT (after options 12(e) and 12(f) 

FI 0 . 0 , 
15 

14 

KODEC, UPC. INTC , HOVE (use with CARD 12 and 

3 OX , 


and 13) 

415 

14 A 

HPCII , t;SEC, IB, LB (if NCDEC<0 on CARD 14) 

415 
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Card Ho. 

Variable 

F orna t 

HL 

XPIV, YPI V , ZPIV, nx, II Y , HZ , ROT (if IJOVE=l 
on CARD 14) 

7 F 1 0 , 0 

15 

TIIETA2, TIIETA1 (only if NODES <0 on CARD 11) 

2F10 .0 

16 

NPC , INTC , KURV, NPTIP, NODES, KPS, NTS 
(special tip patch) (only if I1AKEXD 
on CARD 10) 

3 5X , 
315, 
10X, 
315 

16 A 

(S(I), Y(I), Z(I), 1 = 1, HPT IP (only if EURV > 1 
on CARD 16) 

3F10.0 





Card No. 

Var 1 abl e s 

Foruat 

17 

X (wake grid plane stations) 

F10 .0 

18 

NODE, NPC, INTC, HARE 

3 0 X , 
415 

19 

I DEM TV/, I FLEX'..', I DEFY/ , V/NAHE (wake card) 

315 , 5X 
6A4 

20 

KV/PACII , KW SIDE, KWLINE, KY/PAH1, EY/PAN2 , 
INPUT, NO DEV S , IDV/C, IFLXL, DTKET 

915 , 
F10.0 

21(a) 

SV/PY, SV/PZ, AX (if INPUT=1 on CARD 20) 

3F10.0 

21(b) 

SWFX, SV/PZ , AY (if INPUT* 2 ON CARD 20) 


21(c) 

SWPX, SY/PY, AZ (if INPUT=3 on CARD 20) 


21(d) 

SV/PX, SWPY, SVPZ (if INPUT** 4 on CARD 20) 


22 

NO DEW C , NPC, INTC 

3 OX , 
315 

23 

VIN, VOUT ... (if IDENT¥=4 on CARD 19) 

SF10.0 
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24 

25 


Surface S t reamlin e Inp u t S umm ary 
Variables 


£oriiat 


F, KP , NS (compulsory input if IBLTYP=1) F10.4. 

(Place one card no. 24 for each streamline) 215 

F , KP, NS (end of surface streamline data) 2F10.4, 

215 



B oundary La y er In p ut Summar y 
Variables 

PsNB, TRIPUP, TRIPOP, XPRINT. XSEIP 

(CARD 26 only present if NVPI>0 on CARD 4(a)) 


Fo rm at 

5F10.0 


Off-Body Velocity Scan I npu t Summ ary 

Card No . Vari ab le s 

27 HOLD, I2EET, NEAR, INCPRI, INCPRJ, INCPRK 

(Start of each scan box. Finish the set with 
a blank card) 


F or ma t 
615 


23 

XO, YO, 

ZO. HP (if M0LD=1 

on 

CARD 27) 



3F10.0, 

15 

29 

XI. Yl, 

Zl, NP1 (if HP > 1 

on 

CARD 28) 




2 9 A 

(ALI(I) 

, 1=1, 1 NP1 | ) (only 

if 

IIP 1 < 0 on 

CARD 

29) 

8F10 .0 

30 

X2, Y2 , 

Z2 , HP 2 (if NP > 2 

on 

CARD 28) 



3F1 0 . 0 , 
15 

3 0A 

( AL2 ( I ) 

, 1=1, I NP2 | ) ( on ly 

if 

N P 2 < 0 on 

CARD 

30) 

8F10.0, 

31 

X3, Y3 , 

Z3 , NP3 (if HP=4 

on 

CARD 28) 



3F10 . 0 , 


3 1 A ( AL3 ( I ) , 1-1, t NP3 I ) (only if KP3<0 on CARD 31) 


15 

6F10.0 
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£££iL_Ii£ • Variables 

3 2 XI, Yl, 21, R01, HU , THETAl , T1IETA2 

(if HOLD* 2 on CARD 27) 

33 X2 , Y2 , Z2 , R02, RI2 (if MOLD=2 on CARD 27) 

3 4 NAL, HTIIETA, NR AD (if IIOLD=2 on CARD 27) 

3 4 A (ALl(I) , 1 = 1, | NAL I ) ( if NAL< 0 on CARD 34) 

34B ( ALTIiiET ( I ) , 1*1, | NTH ETA I ) ( if NTH ETA <0 on 

CARD 34) 

3 4C (ALRAD(I), 1 = 1, I NRAD I ) ( if NRAD< 0 on CARD 34) 


E£isii 
7F10 .0 

5F10.0 

315 

GF10 .0 
GF10 . 0 


SF1 O.C 




Off- 

•Eody 

Streamline Input Summary 


Card No. 




V a r i a b 1 e s 

F o rua t 

35 

USX, 

RSY, 

RSZ, 

SU, SD, DELS, NEAR 

6F10.0 


(one 

card) 

card 

per 

streamline; finish with a blank 

15 




7 , 3 He ic r ipt j. o n of Inp u t V a riable s 

Throughout the description below, the following notes apply, 
(i) All integers are right adjusted. 

(ii) Limits on the value of a variable and default 
values are quoted for two versions of the VS AERO 
program: VSAEnO-lOOO^VSAERO-3000 where these 

differ. 

7.3.1 Dasic Inp ut 
£ASH_ii Case Description 

Column s YiLr.i.ab_l e_ Desc r iption Foruat 

1-80 TEXT Alphanumeric text identifying the 20A4 

case 


CARD 2: Control Card 


Col mans 


1-5 

0 Prints all input data 
except the PATCH GEOMETRY 

1 Prints all input data except 

the detail coordinates of the 
PATCH GEOMETRY input. (This 
option is useful for cases with 
large input decks: it still 

displays the information on 
NODE CARDS, SECTION CARDS, etc.) 

2 Prints all input data. (Warning: 

a complicated case may give a large 
amount of printout from the PATCH 
GEOMETRY. For such cases remember 
to request a reasonable output line 
limit; e.g., LGO (PI,=7777 ), or use 
I PR 1=0 or 1 as an alternative) 


Va r i. ab l.£ Va l_ue Fo rm a t 

II'-RI Input data print control 615 
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CARD„2i. C ont inue d 

Va r i.£b_l£ Va lu e 
6-10 IPRLEV 


0 


1 


Des c ripti on Format. 1 

Output print control 

Basic print level (see 
8.0) All other options 
below are additional prints 

Panel corner points arc 
printed (IPRG0IK1); also, 
the wake print is activated 
at the level IPK.VAK-1 (see 
CARD 2 A 


2 


3 


4 


11-15 IPRESS 


0 

1 


N 


The doublet solution, 

|i^( 4nV 00 P.EFL) , is printed 
(IPRSOL-1). Also, the 
neighbor information is 
printed at the level IFRHAC-1 
(see CARD 2 A) 

The P-'Poo values are added to 
the basic printout of panel 
velocities and pressures 
(IPEPPI=1, see under CARD 2A) 

Corner point analysis results 
are printed at the level 
IPilCPV=l (see CARD 2A) 

Gives the option to read in 
the print control variables 
IPRGOi:, IPRNAB, IPRWAK, 

IPRCPV, IPRPPI. Requires 
CARD 2 A to follow CARD 2 

Controls frequency of printout 
in the wake iteration or time- 
stepping loop 

Prints only at the last step 

Prints at every step (avoid 
this if using a large number 
of time steps or wake relaxa- 
tions) 

Prints at every Nth step (includ- 
ing the first or last steps) ac- 
cording to the IPRLEV option 
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CASl£_2i Continued 

Columns V ar i able Valu e 



Description 


Format 


16-20 IISTOP 0 Complete run through the 

code 


1 Calculations stop after 

GEOllIN subroutine. This 
sets 1 PRGOU ,b 1. (The input 
file only requires basic 
data and patch geometry for 
cases with I.1ST0P=1 or 2) 

2 Calculations stop after the 

SURPAN subroutine. This sets 
IPRCO!l=2 . The panel geometry 
file is formed (TAPE 7) and 
can be saved for plotting pur- 
poses 

3 Calculations stop after the 

V/AKPAN subroutine. This sets 
IPRHA3-1 and IPRWAK=2. The 
matrix of influence coeffi- 
cients is not formed if 
NSTOP=3 . The geometry file 
(TAPE 7) now includes the 
initial wal:e geometry* (The 
input file must now include 
the WAKE INPUT) 

4 This is a complete run through 

the code but a RESTART file is 
formed after subroutine ANAL I Z 


Note : The print control values set on IPRGOI1, IPEHAF and IPRWAK 

according to t ii C STOP values above, can be overwritten 
using the IPRLEV options on this card (e.g.,see CARD 2A). 


21-25 HSTART 


Restart control (see 6.3) 


0 Regular first run of a case 

3 Program restart for further 

solutions and wake shape 
iterations and/or boundary 
layer calculations (IELTYP=0) 


7 C 



£A&P_ 2: . 

Colu mns 

21-25 


26-30 


CARP_lAi 

Columns. 

1-5 


Coat lnned 

X*£i»kie iaine P«*££i£t ion Format 


N START 4 


5 


MODIFY 


0 


>0 


2 


Program restart for more sur- 
face streamlines aad bouadary 
layer calculations (IBLTYP=1) 
Future option 

Program restart maialy for 
off-body velocity scans aad 
off-body streamliaes usiag 
the earlier solution 

Controls details of a re- 
start run (i.e., MSTART>0) 

No change in the basic condi- 
tions except on CARDS 1 
through 4 


All the BASIC INPUT must be 
repeated as for the origiual 
run but with the required 
changes for the restart case 

The WAKE INPUT must be included 
with the required changes 


Additional Print Control 

XasUkU X i. las. 


(Only if IPRLEV=5 on 


CARD 2) 
format 


IPRGOM 


0 


Controls printout of 415 

panel geometry (see 8.0) 

Print off 


1 Panel corner points printed 

for all panels 

200+N Panel corner points printed 

for panels on patch N only 

2 Panel control points and 

unit normal vectors printed 
for all panels 

400+N As for IPRGOM-2 but for panels 

on patch N only 
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Card 2A 

Continued 




Columns 


Value 

Description 

Format 

1-5 

IPRGOII 

-i 

Prints basic points (input 
and generated) for defined 
sections on all patches. 



- ( 2Q0+N) 

As for IPRGOII*- 1 but 
patch N only 

for 


Note ; Negative IPRGOII values are intended for use in runs with 

'NS TCP* 1 (see under CARD 2) for the purpose of 
identifying strings of basic points for copying into 
later defined sections. 


6-10 IPRNAD 

0 

1 

2 

11-15 I PRY/ AH 

0 

1 

2 

3 

16-20 IPLCPV 

0 


Controls printout of panel 
neighbor information (see 3.0) 

Off 

Inform ation is printed for 
panels on patch edges, panels 
at wake-shedding lines and 
panels which have failed to 
find one or more of their 
neighbors 

Prints neighbor information 
for all panels 

Controls printout of wake 
data (see 3.0) 

Off 

Prints wake-shedding infor- 
mation for each wake column 

As for IPIIY7AK = 1 plus details 
of wake line geometry 

As for IPRV/Ak= 2 plus wake 
panel doublet values 

Controls printout of panel 
corner point analysis 

Off 


£0 



CARD 2 A : Continued 


Columns 

16-20 


21-25 


CARD 3 : 
Coluuns 
1-5 


6-10 

11-15 


16-20 


¥a.r_i a_b Value 

IPRCPV 1 


2 


Description Forma t 

Prints x, y, z, vx, vy, vz 
V and CP for each panel 
corner point 

As for 1 plus panel corner 
point doublet and source 
values 


IPIiPPI Controls printout of P-'Poo 

values with the velocity 
and pressure data at panel 
centers 


0 Off 


1 On 


HOPE CARD 


Variable Value 


De s c r ip t i on 


F or na t 


LODE 


1 Steady calculation. 815 

Uses option (a) on CARD 4 


2 Unsteady' calculation. 

Uses option (b) on CARD 4 
and requires CARD 6 A 

N PM MAX Upper limit on number of 

panels the user expects 
to be generated (default 
1000^3000) 


NRDMAX Limit on the block size in 

the blocked Gauss-Scidel 
solution iteration. (Default 
value — also the upper limit-— 
is 14CW90) 


-N A negative number allows the 

user to specify the block 
sizes in I! blocks. Requires 
CARD 3 A to follow 

ITGSL.X Limit on the number of Gauss- 

Seidel iterations for a solu- 
tion (default = 20) 
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CARD 3 : Coat inued 

C olum ns Var i a b 1 e Value Descriptio n Format 


21-2S I HEDGE 

26-3 0 N DUE 


31-35 MSPHAX 


36-40 NPCNAX 


0 Vortex line merging 

future option 

Number of subpanel inter- 
vals used on a near-field 
wake panel whose IDENTU > 1 
during the evaluation of 
wake influences in PIIIKAT 
subroutine. These sub- 
panels are used only in 
the streauwise direction 
here. (Default = 10) 

Limit on the number of sub- 
panels per panel used on 
near-field panels on a wake 
in the velocity (VEL) routine. 
(Default = 25) 

Limit on the number of 
Predictor -'Corrector cycles 
in the steady (IiODE=l) wake 
relaxation in subroutine 
u AIIEEL . (Default = 2) 


Note : NED II AX and ITGSI'.X are only active if the number of 

panels exceeds 320-'635 . For a stialler number of 
panels, a direct solver is used. 


.CARD 3 A : User-Specif i ed Clock Sizes in Bl o cked Gauss-S e idel 

Routine ( DUBSOL) (Cn ly if NRBIIAX negative on CARD 3) 


C olumns Var .ijtb l.e. ya lu e 


De sc ripti on Fo rm at 


1-80 NROWB(I), 

1 = 1. N 


V here N=|NRBIIAX| on 1615 

CARD 3. Number of rows 
in each block of the 
matrix 
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CAR D 4 : N uub er o f Step s . Use one of two options depending on 

the value of MODE on CARD 3. 


0PTI011_ia.li ( II0DE=1 ) 

Columns Variab le Value 

1-5 NUIT 0 


N 


6-10 NVPI 0 


Descrip t io n poriaa t 

One pass through poten- 
tial flow calculation; 
i.e., rigid wake 

Number of wake shape 315 

iterations per viscous^ 
potential flow iteration 

Potential flow solution 
on 1 y 


>0 Number of vi scous^po t en t ia 1 

flow iterations. CARD 26 
must be included in the 
input deck 

11-15 IBLTYP Type of boundary layer 

procedure if NVPI>0 


0 Stripwise infinite swept 
wing. Requires CARD SET 4A 
to input the required 
strips for this analysis 

1 Two-dimensional procedure 
along calculated streamlines. 
SURFACE STREAMLINE INPUT 
(CARD SET 24) must be used 
for this analysis 


OPTION (b ): (MODE= 2l 


Columns 

Variable 

Value 


Description 

1-5 

NT 


Number 

of time steps 

6-10 

NIIC 

23 

Number 

of half cycles 


F orm a t 
215 
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ORIGINAL 
OF POOR 



CARD SET 4A . Specifies patch columns for stripwise boundary 

' layer calculations. Only if NVPI>0 and IELTYP=0 

on CARD 4(a) . (IlODE=l) . 


£AR D._4A ( i ) : Number of Sets of Columns 

Columns Variabl e Value Description Format 

1-5 NPSETS I£NPSETS£5 Number of sets of patch 15 

columns 


CARD_4Aiiil 

Columns. Va r ia^b _le^ VaJLji®. D_e^£r t. J.p^n F orma t 

1-5 NPCIIBL Patch number (note: the 1615 

patch should have I DENT® 1 
on CARD 10) 

6-10 NDCOL l£NECOL£14 Number of columns on 

this patch on which the 
stripwise boundary layer 
calculation is required. 

At least two columns should 
be specified (preferably 
the first and last). If 
KECOL® 1 , then all the 
columns on this patch will 
be analysed 

11-80 ( IIO L ( I ) , 1 = 1 , Column numbers (numbered 

NBCOL) relative to the start of 

this patch) only if NEC0L>1 


Note : Use one card 4A(ii) for each set. If more than 14 

columns are required to be analysed on one patch, then 
an additional set can be used with the same patch number. 

If NSTARDO and MOBIFY=0, then proceed to end of BASIC 
INPUT. 
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CAItD_Jij S ymme try C ard 

Va riable Valu e 
1-10 RSYli 0.0 

1.0 

11-20 I! GPU 0.0 

1.0 

21-30 RNF 

31-40 RFF 

41-5 0 RCOUE 

51-60 SCLP.ES 

61-70 TOL 


De s cription 

Symmetrical case (aero- 
dynamic and geometric 
symmetry about the j =0 
plane ) 

Asymmetrical case 

Free-air case (no 
ground plane) 

Ground plane present on 
the z = 0 plane 

Radius of near-field 
factor; default = 2.5. 
(Cased on the panel mean 
v/ i d t h ) 

Radius of far-field 
factor; default = 5.0. 
(Cased on the panel 
mean width) 

Core radius on vortex 
filaments (based on CEAE-' 
2.0); default = .05 

Gauss-Seidel Residual 
Limit as a percentage of 
the maximum doublet size 
in the solution; default 
= .2 

Tolerance limit for the 
test of proximity of a 
point to a panel edge. 
(Factor on the panel side 
length); default = .2 


Format 

7F10.0 
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CARD 6 j__ 

Onset Conditions 


Coi umn.s 

Variable 

Value 

Description Format 

1-10 

ALDEG 


Incidence of x-axis in 4F10.G 

degrees 

11-20 

YAW DEG 


Yaw of x-aiis in degrees 

21-30 

REACH 


Each number 

31-40 

VMOD 


Onset flow velocity magni- 
tude (default = 1.0 if 
VflOD < 10"° 

41-50 

COHFAC 


Compressibility algori t hn 
factor (default = 0.0). 
Second-order correction 
being checked out 


£AAD_6 Aj (Only Pr e sent if UODE =2) 


unn s. 

Variable 

Value 

Description 

1-10 

ALEAR 


Asp 1 i t u d e of motion 
(degrees) 

11-20 

RFREU 


Reduced frequency, 
nfCEAR^V 

OO 

21-30 

I£X \ 
UY l 
iiZ ) 
\ 


Pivot axis unit vector 
tk rough reference monent 
point (default 0.0, 1.0, 

0.0) 


Format 
5F10. 0 


CAR D 7 ; R efe r en c e C on dition 

C olumns Variable Va .1 u£ Description Form at 


1-10 CBAR Reference chord used 6F10.0 

for normalizing pitching 
moment, (Note: reference 

length ( REFL ) used inside 
the code for the unsteady 
node (!IODE=2) and also 
for normalizing the geo- 
metry is CDAR^2.0) 
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CAI1D,.„1i, Con tinue d 

Col unn s YiLI.i.Jlb Jjb Val ue 

11-20 SREF 

21-30 SSPAK 


31-40 

41-50 

51-00 


RLPX ) 
niiPY \ 
ShPZ ) 


D es cr iptio n 

Reference area 

Senispan used for 
normalizing rolling 
and yawing moments 

Coordinates of the 
reference moment point. 
If I.I0DE=2, this is a 
point on the pivot axis 


Format 


CARD_Gj S pecial Opt ions 

C o_l unn_s Va.r i. a.b Xs. Value 

1-5 IIOR SET 


6-10 NVORT 


11-15 NPASUI2 


16-20 JETPAH 


De scri p tio n For m at 

Allows non-zero normal 515 

velocities to be set. 

Requires CARD SET 8A if 
NOR SET) 0 . HORSET = number 
of sets of panels (limit 
= 200^300) 

Allows a free vortex to 
be specified. NVORT = 
number of segments on 
vortex. Requires CARD 
SET 8D if NVORT >0 (limit 
= 20). This option is 
not complete at this time 

Allows sets of panels to be 
identified for a separate 
accumulation of force and 
moment data. Requires CARD 
SET 8C if NPASUfl > 0 . NPASUM 
= number of panels sets 
(limit = 200^300) 

Allows panels to be identi- 
fied inside jet regions of 
high or low energy to make 
correction for total head 
difference in the analysis stage. 
Requires CARD SET 8D if JETPAN>0. 
JETPAN = number of panel sets 
limit = 200^300) 
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GAUD 8: Continued 


Columns. 

Variable 

Value 


Format 

21-25 

NBCUGE 


Allows the automatic 



panel neighbor infor- 
mation to be overwritten. 
Requires CARD SET 8E. 

(Note: an initial run 

with I»ST0P=3 and with 
IPRLEV=2 on CARD 2 
gives a quick look at the 
automatic panel neighbor 
information. The present 
option allows that data to 
be changed on a panel by 
panel basis) 

Note : The action of NOR SET and JETPAN also affects the panel 

neighbor relationships so that doublet gradients are not 
attempted across boundaries of different regions. 


CA RD SE T 8 A : Only present if N0 RSET> 0 on C ARD S 


Columns 

Variable 

Value 

Description 

Fornat 

1-5 

NORPCII 


Number of patch con- 
taining the required 
panel set 

515 

2F10.0 

6-10 

11-15 

NGRF \ 
NQRL J 


First and last rows of 
panels, inclusive, cover- 
ing the required set of 
panels, see Figure 31 
(default gives full set 
of rows on this patch) 


16-20 

21-25 

HOCF \ 
NO CL J 


First and last column of 
panels, inclusive. (De- 

fault gives full set of 
columns on this patch) 


26-35 

VNORIi 


Specified normal velocity 
for the set of panels identi- 
fied above. Positive out- 

wards from the surface 





FIRST AND LAST 
ROWS IN THE SET 

\ f 


FIRST AND UST 
COLUMNS IN THE SET 


SIDE 3 


IDE 1 „ 


PANELS WITH SPECIFIED NORMAL 
VELOCITY, VNORM, ON CARD 8A 
FOR PATCH NORPCH, VNORM IS 
CONSTANT OVER THIS PANEL SET, 
(SIMILAR TREATMENT FOR CARDS 
8C AND 8D) 


SIDE 2 


Figure 31. Identifying a Set of Panels on a Patch, 









Of POOR QUALITY 


CAuD_8Ai Continue d 

C ol umn Var i abl e Va..lu£ Descr iption For nat 

36-45 ADUC Applied doublet value-- 

allows an initial doublet 
solution to be applied at 
the outset 


No te ; Use one card per set of panels. Total number of sets 
= HOE SET on CARD S. 


C ARD SET SD : Only Present if NVQRT > 0 on CARD 8 . 

CARD_SDiil 


Col umns 

Var i ab 1 e 

Value 


Forma t 

1-10 

VORT 


Vortex Strength 

F10.0 


CARD 8D( i i ) 



Columns 

Variable 

Va 1 ue 

Description Format 

1-10 

RXV ( I ) 


Points describing 3F10.0 

vortex line 

11-20 

EYV ( I ) 



21-30 

RZV(I) 




1 = 1. 

NVORT+l ; 

One card per point. 
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ORIGINAL PAGE 11 
OF POOR QUALITY 


CARD SET 

SC : Only 

Present if 

NPASUII >0 on CARD 8. 


Columns 

Variable 

Value 

Description 

Format 

1-5 

NPS pen 


Humber of patch contain- 
ing set of panels 

515 

6-15 

11-15 

NPSRF \ 
NPSP.L ( 


First and last rows of 
panels, inclusive, cover- 



ing the required set of 
panels, see Figure 31. 
(Default gives full set of 
rows on this patch) 

16-2C KPSCF I First and last colanns of 

21-25 NPSCL f panels, inclusive, covering 

the required set of panels. 
(Default gives full set of 
columns of this patch) 

No t e : One card per set. Number of sets - NPASUIi. 


CARD_SET_CDl QnlZ_£resent_if_JETPAN>0_ on C ARD S. 


Columns Variable Value 


Descripti on 


Format 


1-5 

JETPCH 

6-10 

11-15 

JETRF 

JETRL 



16-20 

21-25 

JETCF 

JETCL 

26-35 

36-45 

VIN 1 
VOUT J 


Number of patch contain- 515, 

ing the required panel 2F10.0 

First and last rows of 
panels, inclusive, cover- 
ing the required set of 
panels, see Figure 31. 

(Default gives full set of 
rows on this patch) 

First and last columns of 
panels, inclusive, cover- 
ing the required set of 
panels. (Default gives full 
set of columns on this patch) 

Tangential velocity components 
in the 'jet-wise' direction on 
the inside and outside surfaces, 
respectively, on the jet sheet 
enclosing region. (Normal vec- 
tor on jet sheet points outside) 


91 



CARD SET SIS,: 


Qftly_Pt e sent if NpCIiGE>O wi .pn m CARD, 8 . 


Col unn s 

Variable 

Value 


1-5 

6-10 

KPAN 1 
KSIDE f 


Panel number and the 415 

side of that panel 
requiring a modified 
neighbor 

11-25 

16-20 

NEVHAD l 
NEW SID ( 


New neighbor and the side 
of that neighbor adjacent 
to KSIDE of KPAN, see 
Figure 24 


ii£i£ : One card per set; number of sets = NDCIiGE . 

If NEVMAB = 0 then the program sets NKWSID = -KSIDE. 
Also, it treats the original neighbor in the same 
way. Thus, neighbor relationship is cut completely 
across side KSIDE, 


7.3.2 £ ATCII^GEO IIETR Y^INPUT 

Repeat the following cards for each coaponent^pat cli as 
appropriate . 


£ ARD_9 J £o^R£££££_£jL£d • 


Columns Variable Value 


Description 


Format 


1-10 

CTX) 

Location of component 

5F10.0 

11-20 

CTY S 

coordinate system origin 


21-30 

CTZ ) 

in the global coordinate 
system, see Figure 7 n 


31-40 

SCAL 

Scale factor ( I SCAL 1 ) 



applied to component geo- 
metry (default - 1.0). 

If set negative, include 
CARD 9 A 

41-50 TIIET Rotation angle (degrees) 

for the complete component 
(default axis = component 
y-axis) , For an arbitrary 
axis use negative SCAL and 
use CARD 9A 

£i£fJL : Limit on number of components = 10. 




CA RD 9 A: (Only Pr ese nt if SC AL< 0 on CAR D 9) . 

£2.iu.£ili.£ YiL£i.JLlL.l.e Ya.l_ u£ Fo rmat 


1-10 CPX 
11-20 CPY 
21-30 CP2 

31-40 C1IX 
41-50 CIIY 
51-60 CHS 


I 


Two points on the arbi- 6F10.0 
trary axis about which 
the component is to be 
rotated 

(CPX, CPY, CPZ becomes a 
point on the axis and CIIX, 

CIIY, CIIZ is changed to the 
unit vector along the axis 
of rotation) 


CARD 1 0: Patch Card , 

Col uan s Va^_i a.b .le. Va„lue_ De s c r i u t i o n Form at 

1-5 IDliNT 1 Wing- type patch, ana ly- 41 5 , 6 A4 

sis includes sectional 
quantities. Allows 
IELTYP=0 boundary layer 
analysis (see CARD 4(a)) 

2 Body- type patch, analy- 
sis excludes sectional 
quantities. Only IBLTYP=1 
boundary layer analysis 
allowed on this patch (see 
CARD 4(a)) 

3 Patch with Neumann boundary 
condition panels (single 

s heet ) . Mo boundary layer 
analysis allowed on these 
patches at this time 

6-10 MARE Switch for automatic patch 

generator 

0 Regular patch input. Requires 

CARD 11 and CARDS 12, 13, ... 

w he re applicable 

+Ii Automatic patch closing side 3 

of folded patch, 14; see Figure 
22 

-II Automatic patch closing side 1 

of folded patch, K (Continued) 
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CA2D_.ji.0i., Con t inue d 

N ot e: (i) ilAKE\0 requires CARD SET 16 only to follow CAR!' 10. 

(ii) If 11AEE\0, the upper and lower surfaces of the 

folded patch, | MAKE I , should have the sane number 
of panels. 

Columns Va ria ble Va._l.ue J2£.A££i£jtion Fo rua t. 

11-15 ROMP COMPONENT number to which 

this patch is assigned (C 
or blank automatically 
assigns the patch to the 
current component number). 

Note: the f i r s t patch 

on a new component should 
not be assigned to a dif- 
ferent component. See 4.4 

16-20 KLAS5 Assigns an ASSEMBLY number 

to this patch — these can 
be in any order (default 
i s 1 ) 

Note.: Patches on different COMPONENTS can be contiguous 

(but watch the component transformations on the CARD 9 
information); whereas patches on different ADSENFLir S 
are not allowed to connect in the panel neighbor infor- 
mation. 


21-80 PLANE Text for Patch Identi- 

fication (optional) 

No te: (Limit on Number of Patches = ICO; 

Limit on Panel s-'Pa t ch = 5 00^1000; 

Limit on IIROW or NCCL for a patch = 100.) 


CAP.D 1 1 : Section Card 


Columns 

Variable 

Value 

Description 

Format 

1-10 

STX) 


Location of the 

section 

6F10. 0 

11-20 

STY \ 


origin in the component 

415 

21-30 

STZ ) 


reference frame 
F i g ur e 11) 

(see 


31-40 

SCALE 

20 

Scaling factor 
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ORIGINAL PAGE IS 
OF POOR QUALITY 

CARD 1 1 : Continued. 


Columns 

41-50 

Variable 

ALF 

Value 

Description 

Section pitch angle in 
degree s | 

For ua t, 

i 

5 1 — 0 0 

THETA 


I 

Section orientation / 

angle in degrees | 

^ Sec Figure 

111 

G 1 - 6 5 

IlIuODE 

A 

Type of section input 



0 Copies the previously de- 

fined section as originally 
specified; i . e . , basic co- 
ordinates before scaling, etc. 
If the original section was de- 
fined using INIiOPE=4, then new 
section is located using incre- 
ciental vector ( STX, STY, GTZ) re- 
lative to the original location 


-II 


1 


2 


4 


Copies the basic coordi- 
nates of section K (abso- 
lute subscript) 

Input y , z coordinates of 
SECTION, Requires Option 
(a) on CARD SET 12 


Input x , z c o o r d i u a tea of 
SECTION. Requires 
Option (b) on CARD SET 12 

Input x,y coordinates of 
SECTION. Requires 
Option (c) on CARD SET 12 


See 

Figure 

12 


Input x,y,z coordinates of 
SECTION. Requires Option 
(d) on CARD SET 12 


Generates coordinates on 
a NACA 4-digit section. 
Requires Option (e) on 
C AR D S n T 12, sec Figure 17 


6 Generates coordinates on 

a semi-ellipse section. 
Requires Option (f) on 
CARD SET 12 (See Figure 12) 


7 Generates coordinates on 

a biconvex section. Re- 
quires Option (e) on CARD 
SET 12 (sec Figure 17) 

9 5 



card 1 1 ; 


Coat iniied , 


Q'SLkRI±Ii$ Var AaMe. YiLl u£ 

61-65 imHode 8 

{ cont inucd) 


12 


66-7C NODES C 


1 

2 


3 


A 


N e g a t i v e 
(1 thru 5) 


!2iLiL££A2ii.£2i F o riaajt 

Generates a conpiete 
ellipse. Requires Option 
(f) on CARD SET 12 (see 
Figure 19) 

Polar coordinate input. 

Requires Option (g) on 
C AR D SET 12 (see Figure 
13) 

First or interior section 
on a patch 

End of a spanwise region 
within a patch with con- 
tinuous (1) or discon- 
tinuous (2) slope on the 
spanwise generators onto 
the nest spanwise region 

This section completes the 
present PATCH 

This section completes the 
last patch on the present 
COMPONENT, (The next patch 
will start a new component) 

This section completes the 
last patch in the configura- 
tion. This must be present 
on the last section to ter- 
minate the PATCH GEOMETRY 
INPUT 

If MODES is set to a negative 
value, the present section is 
rotated about the S.C.S. 
x-axis to form a part or a 
conpiete body of revolution. 
Requires CAEd 15 to follow the 
section coordinates (see 
Figure 21) 



CA R D_1 1 : C o a t. .i n u e d , 

Va riable Value 
71-75 KPS 

0 

>0 

76-80 HITS 


0 

1 

2 


Note: Limit on the number of 


2£JL££i£Li. on F or mat 

Only active if NODE5XO 

Iianual paneling in span- 
wise region just completed; 
defined sections are used 
as panel edges in this region 

Number of panels to be 
generated in the s panwise 
region just c oap 1 eted 

(Only active if NCDE5X0 
and NPS > 0 

Form of spanwise interval 
spacing for the generated 
panels in the spanwise 
region just completed. 

See Figure 15 

Full cosine spacing with 
s .m.fi Her panels near t he 
beginning and end of the 
re g i on 

Half-cosine spacing with 
smaller panels near the 
beginning of the region 

Half-cosine spacing with 
smaller panels near the 
end of the region 

Equal spacing throughout 
the region 

defined sections: 500. 


CARD SET 12: Section Definit i on . (Present if IIIIIODE>0 on CARD 

11, see 4.5.5). I r. each set use one card per point. 
Insert NODE CARD ( S ) 14 to control paneling and to 
complete a set. Use one of the options below 
•depending on the value of IHIiODE. Use one card 
per point if I N 1*0 Da — 1 , 2, 3, 4 or 12; see Figures 
12 and 13. 
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CA XD SET 1 2: Continued 


OPTION (a): (IHIK) DE=1) 


Columns. Va r i ab 1 e Valu e 


D escription 


1-10 

LY 

11-30 

LZ 

21-30 

AX 


y,z coordinates of a 
point on the section. 

The x-stations are es- 
sentially constant (0.0); 
however, local deviations 
in x can be placed in Ax 


opTion_ii,ii ( ni:io Di>2) 

Ssinsas Var i a b 1 e Va. l_u e_ De scriptio n 


1-10 

EX 

11-20 

i:z 

21-3 0 

AY 


x-z coordinates of a 
point on the section. 

The y-stations are es- 
sentially constant (0.0); 
bowever, local deviations 
in y can be placed in Ay 


OPTION (c ) : ( IN11 QDE = 3 ). 

C ol mans Variable Va l_ue_ Descripti on 


1-10 

EX 

11- 

-20 

LY 

21- 

-30 

AT, 


x-y coordinates of a 
point on the section. 
The z-stations are 
essentially constant 
(0.0); however, local 
deviations in z can be 
placed in Az 


OP TION (d) : (I HIIODE -4) 


Columns 

Variable 

Value 

Description 

1-10 

TiZ) 


z, y , z coordinates 

11-20 

LY J 


of a point on an arbi 

21-30 

LZ ) 


trary skewed section 


For ua t 

3F10.0 


Format 

3F10.0 


Forua t 

3F10.0 


Forma t 
3F10.0 




OPTIOM (e): INIi ODE=5^NACA 4-D IGI T S ECTION, or I1I!:QDL = 7^CICG1.VE:.: 



SECTION 



o 1 mans 

Variable Value 

Description 

Format 

1-10 

TC 

Thi ckne s s^chord ratio of 

F 1 0 . 0 . 



section 

15 


11-15 INPUT 


1 Generates y,z coordi- 

nates; x = O.G 


2 Generates x,z coordi- 

nates; y = 0.0 


3 Generates x,y coordi- 

nates; z = 0.0 

Note : Tiie coordinates are generated on a chord of 0.0 to 1.0. 

They start at the trailing edge on the lover side and 
finish at the trailing edge on the upper side. (Dec 
Figure 17.) 


CARD SET 13 rust follow this card. 


OPTION ( f ) : IHIiCDE=6 

C o lua he Variable 


1-10 


SEuI-ELLIPSE_or_INE 0 jiE 58 ^_FULL_ELLIPGE. 


Value 

Description 

Fornat 


Height of 'vertical' 

F10.0. 


seaii-axis; see Figures 

15 


10 and IS 



11-15 INPUT 1 


Generates y , z coordi- 
nates; x = 0.0 


2 Generates x , z coordi- 
nates; y = 0.0 

3 Cer.erates x , y coordi- 
nates; z = 0.0 


Note: CARD SET 13 uust follow this card. 
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OPTION (r. 

.Ll— ilHiififlfi 

=12) Polar 

_C o££din a t. e._I np,u t. . 


Col_ ULili S. 

Variable 

Value 

££JL£JLi£iA£ii 

V o r u a t 

H 

1 

O 

EX 


x- .station 

3110. 0 

11-20 

RAD 


Radius vector 


21-30 

TEXT 


Orientation (degrees), 
see Figure 13 


CARD SET 

.13: (Used 

After Option 12(e) or 12(f) . 


Col unn s 

Variable 

Value 

£££££A£.JL£££ 

Fornat 

1-10 

XRD 


Chordwise station in the 
range 0.0 to 1,0 at the 

FI 0 . 0 , 
15 


end of a cliordwi.se region 
on a generated section. 

Place a negative sign on 
XRR if tiie region end occurs 
on the section underside. 

See Figures 17, 18 and 19. 

on the last CARP 13 of a 
set, XR1';= 1.0 for a c o r.i p 1 e t e 
section except in the case 
when INfiODE=6 : in this case, 

the last region ends with 
XP.D=0 . 0 

11-15 HINT Number of basic point inter- 

vals to be generated in the 
chordwise region just con- 
pie ted. Default value is 
70 

£i£l£ : A NODE CARD (CARD 14) must be placed after each CARD 13. 

The TERMINAL NODE CARD (with NC'DEC = 3 ) nust be placed after 
the last CAP.D 13 in the set. 


CASD_14i Chordwise Node Card . 

Columns. Va_r i.a b 1.0 YxAue De sc r ip tion 

30-35 N0D2C ll Terminates a chordwise 

2J region having continuous 

(1) or discontinuous (2) 
surface slope onto the 
next chordwise region 
IOC 


Foroa t 
3 OX, 41 5 



CARD 14; 

Columns 

30-35 


36-40 


41-45 

46-50 


ORIGINAL PAGE If 
OF POOR QUALITY 


Variable Value 

NODEC 3 

Continued 



4 


NPC 0 


>0 


INTC 



HOVE 


0 

1 


Signifies a TERMINAL 
NODE CARD. This must 
be placed after the last 
point on a section 

Negative values indicate a 
copying routine. Use CARD 
14A to define a string of 
points to be copied over to 
form part or all of present 
section. The last copied 
point is the end of a chord- 
wise region with the corres- 
ponding action according to 
the modulus of the NODEC 
value, see 4 . 5 .7 . 1 

As above but the last point 
copied is not at a chordwise 
region end. In this case the 
CARD 14A must be followed 
either with another negative 
NODE CARD or with further 
basic points 

Manual paneling in the chord- 
wise region just completed 
(i.e., basic points correspond 
to panel corse r s ) 

Number of panels to be generated 
in the chordwise region just 
completed 

Form of chordwise interval 
spacing for the generated 
panels. See 4.5.5 and Figure 
15. Also under CARD 11 

(Only used for basic point 
input; i.e., INMODE-1, 2, 3 
4 or 12) 

No action 

Chordwise region just com- 
pleted is to be transformed. 

CARD 14B must follow 
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.BolZ_of Revolution 
on CARD 11.) 


Columns 


Yllfie 


1-10 THETA2 

11-20 THETA1 



a string of basic points to be 
lows CARD 1 4 if NODEC negative.) 


fiei££i£tion Format 

Patch number 415 

Section number within 
that patch (local sub- 
script ) 

First and last basic 
points on that section 
(local subscript) 


ion . (Only if 


fi££££iRii££ £oraa t 

Pivot point coordinates 7F10.0 


Pivot axis (unit vector) 


Rotation angle (degrees) 
(Positive right-hand 
rotation) 


(Only present if NODES negative 


Description Forma t 

Theta interval for body 2F10.0 

of revolution. This can 
be used to build up a body 
of revolution with varying 
panel density. See Figure 21 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


(LAR D-M X -Ml — &LL&1&1-JEA P Pitch. See 4. 5. 7. 2 and Figure 22. 
(Only present if MAKE\o on card 10) 

£ &l ia& l . XJLLLiMS. IaIr* Sescr.ig.tion Format 


36-40 NPC 


Number of panel 'chord- 35X, 415 

vise' across patch 101,315 


41-45 INTC 

46-50 IURV 


51-55 NPTIP 


66-70 NODES 


Form of panel spacing 
'ehordvise' (0, 1, 2 or 
3 options as under CA&D 11) 

0 Flat patch 

1 Semicircular sections 

2 Semi-elliptical sections 

3 Triangular sections 

Number of points defining 
the tip-edge contour. 
Required only for KURV>1, 
then use CARD 16A 

3 More patches to follov 
this one 

4 This is the last patch on 

a component 

5 This is the last patch in 

the PATCH GEOMETRY input 


71-75 NPS ) 0 Same options as given on 
76-80 INTSJ 0 CARD 11 but not usually 

needed here (the default 
paneling matches that of 
the basic patch) 


CARD SET 16A: Tin Planform Shape . (Only if NPTIP >0 on CARD 16) 


Columns Variable. Xalue 


Description Format 


11-20 S ) 

11-20 Y } 

21-30 Z ) 


Coordinates describing 3F10.0 

the tip-edge shape in a 
local system. See Figure 
22(d). z-0.0 if KURV-2 


Note : Use one card per point. Number of cards 


NPTIP. 
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7 . 3.3 £a fee m IS£U t 


ORIGINAL PAGE IS 
OF POOR QUALITY 


CM2_1 2i__VakerC rid -Elan e . 
CfiiSSUli. Yari.£bl.e Valu e 
1-10 a 


Pe scr i p tion F orna t 

x-station of eacL wake- F10.0 

grid-plane; one value 
per card. First station 
oust be upstrean of all 
separation points 


The set of x-stations must 
be terminated with a NODE 
CARD (CARD 1C) with NCPE=3 . 
Intermediate NODE cards 
(with NCDE=1) may be inserted 
to generate intermediate sta- 
tions as described below. 

(See also Figure 25) 

Note: Limit on Number of Uake-Gr id-Plane s * 31. 


CAR £_ 1 2 i NOD E C ARD . 

Columns Variable Value 

30-35 NODE 1 

3 

36-40 NPC 0 

>0 


Description Foraa t 

Intermediate node card 30X.4I5 

may be placed after any 
CARD 17 except the first 
and the last 

Teruinal node card placed 
after the last CARD 17 

Manual intervals: the set 

of x-stations just completed 
is used di recti y as a set of 
wake-grid-plane locations 

Automatic intervals: NPC 

intervals are generated in 
the region between the pre- 
vious two x-stations 
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CARD IS 
C o l mans 
41-45 


4(5-40 


Note : 


; Coat inued . 

Variable Value 


INTC 


0 


1 


2 


3 

HARK 0 



Description 

(Only active is NPOO.) 
Controls tLe spacing of 
the automatic intervals. 
Spacing options are the 
same as described in 
4. 5. 5. 3 and Figure 15 


Format 


Full cosine spacing with 
smaller intervals near the 
beginning and end of the 
r e g i on 


Half-cosine spacing with 
smaller panels near the 
beginning of the region 

lialf-cosine spacing with 
smaller panels near the 
end of the region 

Equal spacing throughout 
the region 

Usual 


1 Use on one x-station 

to mark downstream end of 
region of interest. If 
ilARR is not used the de- 
fault is two stations 
before last one. (Uake 
downstream of this point 
will not be printed or passed 
to plot file and will not 
receive detail relaxation 
calculation. ) 

Repeat the following cards (19 through 23) as appropriate 
for each wake . 
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gAIB-l&L— I&fc*-gig4 

fi &Ias &i . JmjlLmMa 


XaIbjb. 


E&esaI 


6-10 

11-15 

21-80 

Note; 


ORIGINAL PASS IS 
OF POOR QUALITY 


IDENTW 



IFLEXV 



4 


0 


Type of wake 315 , 5X, 

6A4 

Regular wake . Doublet 
distribution is constant 
in strsaawiss direction 


Unsteady wake (continu- 
ously generated and 
1 r a a spor ted a t r e anw ise 
in tine-stepping loop) 


Jet aodel . Doublet distri- 
bution is linear in streua- 
direction with specified 
gradient. Requires CARD 23 


Flexible wake — will be re- 
laxed if wake shape itera- 
tion specified (NWIT>0 on 
CARD 4(a)) 


1 Rigid wake — will reaain 

fixed throughout wake shape 
iteration cycles 


IDEFW 


WNAME 


0 Separation line defined by 

strings of separation panels. 

Requires CARDS 20. 21 and 22 
where appropriate 

1 Separation line defined by set 

of point coordinates on the sur- 
face. (Future option) 

Text for wake identifica- 
tion (optional) 


Liait on Nuaber of Wakes * 10-^50; 

Limit on Nuaber of Wake Coluans * 50^*200; 
Liait on Nuaber of Wake Panels • 1500-^4000 . 
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ORIGINAL PAGE IS 
Op POOR QUALITY 


CARD 20: Separation Line Specification along Panel Edges . 

(See Figure 28) 


Repeat CARD 20 and its CARD SETS 21 and 22 where ap- 
plicable for each patch crossed by the separation line. 


Columns 

Yar iable 

Value 

fiescr i£tj,on 

Format 

1-5 

KU PACE 


Patch number 

715 

6-10 

STvSIDE 

1.2,3 
or 4 

Patch side which is 
parallel to the direc- 
tion of the separation 
line (see 4.7.3 and 
Note (ii) below 


11-15 

El? LINE 

0 

Separation is from the 



patch edge 

>0 Line number (row or 

column) within the 
patch where wake is 
attached 


16-20 KM P Mil 

21-25 in? PAH 2 


26-30 INPUT 0 


First and last wake shed- 
ding panels in the present 
set. The numbering is local 
along the separation line on 
patch KVPACI' . If the sepa- 
ration line extends across 
a complete row or column of 
this patch in the present 
string, then O/PAHl, EWPAII2 
should be set to 0 (default 
option) 

Copies previous wake line 
geometry 


1 Uake line geometry speci- 
fied by y,z coordinates. 
Requires OPTION (a) on 
CARP 21 

2 Vlake line geometry speci- 
fied by x , z coordinates. 
Requires OPTION (b) on 
CARD 21 


3 hake line geometry speci- 

fied by x,y coordinates. 
Requires OPTION (c) on 
CARP 21 
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Continued 


Colum ns 

26-30 

31-35 


36-40 

31-45 

46-55 


YliiaMe 

INPUT 

(continued) 


HO DEI? S 


4 


£®l££i£tion 


Forma t 


Wake line geometry speci- 
fied by x , y , z coordinates 
(global reference). Re- 
quires OPTION (d) on CARD 21 

First or intermediate 
string of panels being 
specif ied 


3 Completes a wake. Only 

INPUT and the following 
variables are active on 
this card. If INPUT) 0. 
then the appropriate 
geometry of the last wake 
line for this wake (CARD 
21) must follow 


5 


IDWC 


IFLXL 


DTIIET 


As for the KODEWS-3 hut 
this completes the last 
WAKE in the input data 

Future option to change 
the type of wake (IDEMTW) 
on a column-by-colucn 
basis 

Future option to change 
the 'flexible' status 
( IFLEXW ) on a line-by- 
line basis 

Option to rotate a wake 
line geometry about the 
local x-axis 


(i) If. in a subsequent case the panel density is 

changed on a patch crossed by the separation line, 
then KWPAhl, KWPAN2 and possibly KWLINE might need 
to be changed also. 

(ii) The 'direction' of the separation line is such that 

the wake-shedding panels 'upstream' of the wake sepa- 
ration are on the left when looking along the line. 
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£AHD_2_lj. S t icaBwisa.Pale-Liae Geometry , (Only if INPUtVo on 

CARD 20) 

Use one of f our options depending on the value INPUT on 
CARD 2 0. OPTIONS (a), (b) and < c ) require local 

coordinates relative to an origin at the separation 
point with coordinates axes parallel to the global co- 
ordinate system. 

OPTION (d) requires coordinates specified in the global 


Descriptio n Forma t 

y,z coordinates of each 3F10.0 

point on a wake line, 
x values are essentially 
constant (0.0); however, 
local deviations in x can 
be placed in Ax 

:es below.) 


OPTION (b ) : II NPUT= 2 1 . 


Columns 

Variable 

l&l&S 

US.lS.Li2 tion 

Forma t 

1-20 

SWPX) 


x,z coordinates of each 
point on a line, y values 

3F10.0- 

11-20 

SWPZ J 


are essentially constant 
(0.0); however, local de- 


21-30 

Ay 


viations in y can be 
placed in Ay 



(Sec notes below.) 


0£XI0N__icli — IINPUT=il . 

Valu e De sc ription Format 

x,y coordinates of each 3 F 10.0 
point on a wake line, 
z values are essentially 
constant (0.0); however, 
local deviations in z can 
be placed in Az 


golSBRS. 

1-10 

11-20 

21-30 


S17PX 1 
SUPYJ 
Az 


coordinate system. 
0PTI0N_ia.Xi (IHP PT-l) . 


Columns 

Variable 

Value 

1-20 

SV'PY ) 


11-20 

srpz ; 


21-30 

Ax 




(See 


(Tee notes below.) 



oPTion^id • 


Columns Variable Injjut 


Description F orma t 


1-10 

swrx 

x , y , z coordinates of 
each point on a wake 

3F10.0 

11-20 

SWPY 

line specified in glo- 
bal coordinates 


21-30 

SVPZ 




Iloies.: (i) Use one card per point for the above options. 

(ii) Exclude the first point on the wake line (i.e., the 
separation point). This is located by the program 
using the information on CARP 20. 

(iii) Wake node cards. CAP.D 22, can be inserted in the 
set of points to generate additional points on 
curved lines, A terminal node card with K0PE'. : 'C = 3 
must be placed after the last point on a wake line. 
See Figure 29. 

(iv) The points input or generated on a wake line are 
not directly related to the wake panels. The 
x-station$ of the wake panel corner points are 
determined by linear interpolation at the wakc- 
gr id-plane s defined on CA '.IT: 17. 


CARD 22: Wake Mode Card . 


Columns Variable Va 1 u e 

30-35 NOLEW C 1 

2 


3 


Description Forma t 

Optional intermediate 3 OX, 3 15 

node card placed within 
a set of wake line points 
in order to generate addi- 
tional points along a curved 
line with continuous (1) or 
discontinuous (2) slope onto 
the next region 

Terminal node card. This 
must be placed at the end 
of a set of points to finish 
a wake line 



GAUD 2 2: 

Columns 

36-40 


41-45 


GAUD 23 

Columns 

1-10 

c 

11-20 
Note : 


Continued . 
Variable Value 

NFC 0 


>1 


IIJTC 


ORIGINAL PAGE if 
OF POOR QUALITY 


fi® s„£r on Forma t 

Ilanual interval 8 . The 
input points are taken 
directly 

Number of intervals to 
be generated in the 
region just completed 

Form of spacing if IJPC > C 
(see Figure 15 ) 


0 Full cosine spacing with 

smaller intervals near the 
beginning of the region 

1 Half-cosine spacing with 

smaller panels near the 
beginning of the region 

2 Half-cosine spacing with 

smaller panels near the 
end of the region 

3 Equal spacing throughout 
the region 


TYP E -4 V a k e Velo city Data. (Only if INDENTI/-4 on 

CAr.r 19 ) 

Var i ab 1 e Value Description For m at 

VI N Tangential velocity on CF10.0 

the underside or inside 
of the vorteu sheet wake 

VOUT As above but on the top 

side or outside surface 


(i) If just one pair of values is given as above, then the 
code will set the same values for all columns on this 
wake. If different values are required on other columns 
of this wake, then I JVC pairs of values must be input on 
the 8F10.0 format continuing onto additional cards if 
necessary. M'.?C is the total number of columns (i.e., 
number of wake-shedding panels) on this wake. 


Ill 




Sfitje: ( i i ) Surface panels immersed in a jet should he identified 

using JETPAt! on CARD 8 and data on CARD SD. This does 
not affect the doublet jt 0.1 ujfc jl o n — only the anal.xs.Ls of 
surface velocities and pressures; it changes panel 

neighbor information and provides the jump in total 
head for C p calculations inside the jet. 



y 


One data card is required for each streamline to be calcula- 
ted and a final card to signify the end of the streamline data. 
Calculation proceeds upstream and downstream from starting point. 


CARD 2< 


Columns Variable l&l&S. JtS.$.SSA2Sl2.E 

1-10 F 0.05iFi0.95 Location of starting F10.0 

point of streamline as 215 

a fraction of the 
panel side containing it 
(measured anticlockwise) . 

A value of F = 0 . 5 v/ill cause 
the starting point to be in 
the middle of the panel side 


11-15 IIP 


16-20 NS 


Panel number of the panel 
through which the stream- 
line is to pass. Caution : 
this must not be a stagna- 
tion pa is el 

1,2.3 Panel side number for 

or 4 streamline entry. 

Normally to be left 
blank. Program will 
select a side . If a 
side is prescribed 
(li.NSi.4), the user must 
ensure that a streamline 
would actually enter the 
panel through that side 
as opposed to exit ina 
through that side 
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O' ,, 

CARP 25 : Last Card of Streamline Input . 


Columns 

l££iabie 

Value 

Description 

Fornat 

1-10 

F 


Must be set equal to 
2.0 if the body is 
attached to a plane of 
symmetry 

F10 .0 
215 

11-15 

KP 


Must be blank (to signal 
end of input) 


15-20 

KS 


Leave blank for normal 



output: set NS-1 to 

shut off all output 
printing (but save the 
plot file; sec 6.2) 
except listing of input 
data 


7.3.5 Bou ndary Layer Inpu t . (Only present if KVPI > 0 

on CARD 4(a).) 


CARD 26 : Reynolds Nunb er, E.t £ . 

££lil£3LS Variab 1 e Value 2®.s.££ilil.io£ Fornat 

1-10 EKE Reynolds nuuber based 5F10.0 


on reference chord and 
free strcau velocity in 
mill ions , VooC^Vx (10“6) 

11-20 TRIPUP Trip location (x-*c). If 

tripping is not desired> 
TRIPUP®! 

21-30 TRIPOP Trip option. TRIPCP=1 : 

this deters the user from 
specifying a trip location 
where the boundary layer 
could not (because of the 
Reynolds nuaber) become 
turbulent. If too early 
a trip location is speci- 
fied. the location is re- 
positioned to corespond to 
the point where R exceeds 
200 e 
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£M2^£i~_ JSfifill&Ml . 



Columns Variable 
31-40 XPRIHT 


41-50 XSXIP 


Va lu e Description Form at 

Loundary layer printout 

0, No crossflow parameters 
will be printed 

1. Crossflow parameters will 
be printed 

Number of integral inter- 
vals to be skipped between 
boundary layer printouts 
(streamline is divided into 
200 intervals* e.g., if 
XSKIP«10, 20 printouts are 
produced) 


7.3.6 Of f o&zJL&l 0 £iii^£a n_I nRUl 


CAR D 27: S can Pox . 

Columns Variab Ijg, Value 
1-5 HOLD 


1 


2 


0 


6-10 iXET 


0 

1 


2 S.2.££A££i££ Forma t 

Shape of velocity scan 615 

volume* see Figure 30 

Skewed bon option: allows 

a single point, points 
along a straight line or 
straight lines within a 
parallelogram or within a 
parallelepiped. Requires 
CARD SET 28. etc. 

Allows points along radial 
lines in a cylindrical 
volume. Requires CARD 
SET 32, etc. 

Stops the scan 

Controls the intersection 
line routine (see 5.1) 

Ac t ive 

Off 
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F orua t 


OF POOR QUALITY 

CARD 27; Continued . 

Co 1 n ans Variab le Va 1 ue De sc ription 


11-15 NEAR 


Controls the near-field 
routine in the velocity 
calculation (see 5.1) 

0 Active 

1 Active only for surface 
panels (wake near-field 
off) 


1 Inactive 


16-20 IHCPRI ) 

21-25 IN CPU J } 

2 6-3 0 INCPSE ) 


Print frequency control 
for places, lines or 
points, respectively 


0 All the velocity scan 

results are printed 

H Prints the results only 

for every Kth plane, 
line or point, respectively 


Kojte.: Use one CARD 27 for each scan box. Finish the set with a 

blank card (i.e., I.iOLD = 0) . 


£ARD_2Sj Fir st C orner of Skewed Box . (Only if KOLD = l on CARD 27) 


Columns 

Variable 

Value 

Description 

gorna t. 

»— 

i 

o 

XO ) 


Coordinates of first 

3F10.0. 

11-20 

YO ) 


corner point of a box 

15 

21-30 

20 ) 



31-35 

HP 

1 

Single point. No fur- 
ther input for this 'box* 



2 Single line. CARD 29 oust 

follow to complete this 
’box' 


3 Lines within a parallelo- 
gram. CARDS 29, 30 must 
follow to complete this 
'box’ 

4 Lines within a parallelepiped 
CARDS 29, 30 and 31 must follow 
to complete this 'box' 
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or 


£4 ££_2.JLL Me ond^fiorn e r „o f , (Only if I'.OLIW 

CARD 27 a n d HP > 1 on CARP 28) 


Columns 

Variable 

Mm 

1-10 

XI) 


11-20 

Y1 > 


21-30 

ii| 


31-35 

IIP! 

n 


D escription 

Coordinates of a second 
corner point of a box 

Number of points (equally 
spaced) along the straight 
line (XO, VO, 20) (XI, Y.1 , 


Forna t 

3F10.0 

15 


21 ) 


— H A negative sign on the nur.»- 

ber of points allows the 
point locations along the 
line to be specified. Ke-^ 
requires CARD 29 A to follow 


CARD 29A; Specified Point Locations alone First Edge of Dor . 

(Only if NP1 negative on CARP 29) 

Co lumn s Variable Value D escription Forma t 

1-8C (AL1(I), C<AL1£1.0 Normalized location of SF10.0 

I*1,|NP1|) each point along first 

edge of box 


CARD 3 0 : Third Corner Poiat of Skewed Box . (Only if HOLE = 1 on 



CARD 27 and 

HP >2 on 

CARD 28) 


Columns 

Variable 

Value 

Description 

Format 

1-20 

11-20 

21-30 

X2 ) 
Y2 > 
Z2 J 


Coordinates of third 
corner point of a box 

3F10.0 

15 

31-35 

NP2 

N 

Humber of points (equally 
spaced) along the straight 
line (XO, YO, ZO) (X2, 

Y2, 7 . 2 )) . These points 
locate the start of each 
scan line within a scan 
plane 




-11 

A negative sign or HP2 
allows the point locations 



along the line to be specified. 
Requires CARD 3 0 A to follow 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


gAR D „lM» Specified Locat ions of Scan Lines . (Only if NP2 



negative 

on CARD 

30) 


Columns 

Variable 


U§.*£.£iR&ion 

Forma t 

1-80 

( AL2 ( I ) , 

0£AL2£1 . 

0 Normalized location 

SF10.0 


1*1, 1 NP2 | ) 


of each point along the 



second edge of box 


CAFD 31 : Fourth Corner Point of Skewed Box . (Only if H0LD=1 

on CARD 27 and KP-4 on CARD 28) 


Columns 

Variable 

Value 

1-10 

X3 


11-20 

Y3 


21-30 

Z3 


31-35 

NP3 

N 


De s c r ip t ion Portia t 

Coordinates of fourth 3F10.0, 

corner point of a box IS 

Number of points (equally 
spaced) along the straight 
line (XO, YO, ZO) (X3, Y3 , 

Z3). These points locate 
the corners of the scan 
places along the third 
edge of the box 


-N A negative sign on HP3 

allows the point locations 
to be specified. Requires 
CARD 3 1 A to follow 


Specified Locations of Scan Planes . 
negative on CARD 31) 


(Only if NP3 


Colum ns Variable Value 


description 


Format 


1-80 


( AL3 ( I ) , 

1 = 1, 1 NP3 | ) 


Q,£AL3£1.0 Normalized location 

of each point along the 
third edge of the box 


8F10.0 
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CARL. 32: 

-First.. Point 
CARD 27) 

for Cylindrical Volume . (If HOLD = 2 

on 

Co l UXuS S 

Ya£i.able 

Value 

Description 

Format 

1-30 

11-20 

21-30 

Y1 > 
El ) 


Coordinates of first 
point locating axis 
of cylindrical volume 

7F10.0 

31-40 

41-50 

UOl\ 

mi J 


Outer and inner radii, 
respectively (see 
Figure 30(b) 


51-60 

61-70 

TXETA1 ) 
THETA 2 J 


First and last azimuthal 
stations (degrees) measured 


from the vertical (see 
Figure 30(b) 


CARP 33: 

Second Point for_ 

CxIlhil£ical_Voluuc . (If !;CLD-2 on 


CARR 27) 


Columns 

Variable Value 

Be s c ri p_ t j on Forma t 

1-10 

X2 1 

Coordinates of second 5F10.0 

11-20 

Y2> 

point locating the axis 

21-30 

22 ) 

of cylindrical volume 
((XI, Yl, Zl) (X2, Y2, Z2 ) ) 

31-40 

L02\ 

Cuter and inner radii. 

41-50 

RI2 f 

respectively, at second 
sta t i on 


£ARD_34i_ 

Point Distribution 

in Cylindrical Volume. 

(If i:old = 2 


on CARD 27) 



Columns 

Variable Value 

2s. IS. £ intion 

Forma t 

1-5 

NAL u 

Humber of stations 

315 


(equally spaced) along 
the cylinder axis, de- 
fining the location of 
scan planes (see Figure 
30(b) 

-K A negative sign allows 

the stations to be speci 
fied. Requires CARD 34A 
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CARD_3 4 : Continued 


ORIGINAL PAGE Si 
OF POOR quality: 


Col man s Va riable Valu e 

C-10 NTLETA H 


-N 


De sc ription Forua t 

Humber of azimuthal 
stations (equally 
spread between TKETA1 
TIIETA2 ) locating scan 
lines within each scan 
plane 

A negative sign on UTHETA 
allows the stations to be 
specified. Requires 
CARD 34B 


11-15 MR AD N 


-N 


Number of points (equally 
spaced) along each radial 
scan line 

A negative sign on NRAD 
allows the stations to be 
specified. Requires CARD 34C 


CAR D-3 4 A : Specified Location of Sc a n Planes in Cyl i ndr ical 

" Volu me (If“oDE=2 _ on CARP 27 and NAL < 0 on CARD 34) 

C olumn s Variable Value 2£S.££AEl.i.££ F orma t 

1-80 (ALl(I) 0£AL1£1.0 Normalized locations 8F1010 

1 = 1, I MALI ) along the cylinder axis. 

(XI, VI, 21) (22, Y2 , 22) 


CARD 3 4D : Sp ecifie d Loc a tion of S c a n Lines in Cylind r ical. 

VoIuL c .1 1 f i:0DE = 2 on CARD 27 and NTIiETA < 0 onCARD 34) 


Columns Variable Value 


Descr i ption 


Form at 


1-80 ( ALTHET ( I ) , 0£ALTI.‘ET 

1=1, IKTEETAI ) £1.0 


Normalized azimuthal 
stations between TI.'ETAl 
and T1IETA2 


8F10.0 
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CAHD.3 4.C i Spe ci fied Point Locati o ns in Cyl indrica.l V o lume . 

(If ;:0DE=2 on CARD 27 and NF.AD< 0 on CARD 34) 


Column s 

Variable 

Value 

Description 

Fo£ua t. 

1-S0 

(ALRAD(I) , 

1 = 1. |NRAD 1 ) 

OiALEAD 
il . 0 

Normalized point lo- 
cations between inner 
and outer radii on each 
radial scan line 

8F10.0 

Note : 

CARDS 3 4A, 34E 

and 34C 

may all be present (in that 

order) . 


CAHD_35i. 

7.3.7 Off- 
Location of 

l£ii_£t£eanliue_ln£ut 

a Starting Point For Streamline Calc 

ul a t ion 

Columns 

Variable 

Value Da scrip t ion 

gorge t 

1 

o 

R5X 

Coordinates of starting 

6F10.0, 

11-2C 

21-30 

113 Y 

nsz 

point 

15 

31-40 

su 

Distances upstream and 


41-50 

SD 

downs treats , respectively, 
for length of streamline 
measured from the starting 
point 


51-60 

DELS 

13a sic length increment for 



steps along the streamline. 

This is a factor on CHAR 
on CARD 7. (The step length 
is adjusted up and down with- 
in the procedure) 

61-65 HEAR Controls the near-field 

routine in the velocity 
cal cul a tor 

0 Active (use if streamline 

passes close to a surface) 

1 Active only for surface 

panels (wake near field off) 

-1 Off (use if streamline is 

well clear of the surface 
and wake ) 

Note : (i) Use one card per off-body streamline. Finish the set 

with a blank card. 

(ii) This procedure is still being checked out. 
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7.4 Input Flow Chart 


7.4.1 Overview 
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7,4.2 Basic Input 


)F POOR 


CARD 1 : TEXT 

zzimz 


CARD 2; IPRI , IPRLEV, IPRESS, MSTOP , MSTART, MODIFY 


IPRLEV' 
=5? 


YES 


CARD 2A: IPRGOM, IPRNAB , 

IPRWAK , IPRCPV, 
IPRPPI 


NO 




CARD 3: MODE, NPNMAX, NRBMAX, ITGSMX, IMERGE, 

NSUB , NSPMAX, NPCMAX 

I 


NRBMAX 
< 0 ? 


YES 


■H CARD 3 A: NROWB(I) 


NO 


CARD 4: (a) NWIT, NVPI , IBLTYP (MODE=l on CARD 3) | 

or (b) NT, NHC (MODE =2 on CARD 3) 


/ IBLTYP=0 \ 
o)N CARD 4(a)?' 




NO 



ICARD 4 A ( i ) : 

NPSETS | 


CARD 4A(ii) : 

NPCHBL , NBCOL, 


(KOL(I) , 1=1, 


NBCOL } 


'MSTART > O' 

& 

MODIFY=0?y 
CARD 2 


YES 


-»-(end OF BASIC input) 


NO 


CARD 5: RSYM, RGPR, RNF , RFF , RCORE , SOLRES , TOL 

* ..." 


CARD 6: ALDEG, YAWDEG, RMACH , VMOD , COMFAC 

n~ 


'MODE=2 Ot 
CARD 3? 


YES 




CARD 6A: ALBAR, REFREQU , HX, 

HY, HZ 

ZZJ 


CARD 7: CBAR , SREF , SSPAN, RMPX, RMPY , RMPZ 

I 


(continued) 
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BASIC INPUT CONTINUED 


CARD 8: NORSET, NVORT, NPASUM, JETPAN , NBCHGE 


/ NORSET> o\ YES 

(ON CARD 8? > » 


/ NVORT > 0 \ 
[ON CARD 8?; 


CARD SET 8A 

NORPCH , NORF, NORL, 
NOCF , NOCL, VNORM, 
ADUB 

\ 

1 


f 



CARD 8B: 

( i) VORT 

CARD 8B: 

: (ii) (RXV , RYV, RZV) 


/NPASUM>0 \ 

Con card 8?, 


CARD SET 8C: NPSPCH, NPSRF , NPSRL 

NPSCF , NPSCL ' 



1 JETPAN>0\ 
ON CARD 8 ?J 


CARD SET 8D: 



JETPCH, JETRF , JETRL, 
JETCF, JETCL , VIN 
VOUT 


' NBCHGE >0 ^ 
ON CARD 8 V 


CARD SET 8E: KPAN, KSIDE , NEWNAB / 

NEWS ID 
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NEW COMPONENT 


7.4.3 Patch Geometry 


OF POOi QUALITY 



324 





DEFINE ANOTHER SET 


7.4.4 Wake Input 


FROM PATCH GEOMETRY 


CARD SET 17/18: WAKE-GRID-PLANES 

SEE FIG. 25 


CARD 19: IDENTW , IFLEXW, IDEFW, WNAME (WAKE CARD) 



CARD 20: KWPACH, KWSIDE, KWLINE , KWPAN1, KWPAN2 

INPUT, NODEWS, IDWC, IFLXL , DTHET 


* INPUT? 
.(CARD 20) 


CARD SET 21/22 



STREAMWISE WAKE LINE COMPLETE 


NODEWS? ' 
(CARD 20) 


IDENTW=4?S 
N CARD 19> 


CARD SET 23: VIN, VOUT 


NODEWS ?\ 
ON CARD 20, 



END OF WAKE INPUT 
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ORIGINAL FAGS IS 
OF POOR QUALITY 


7.4.5 Surface Streamline and Boundary Layer Input 
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6 Of f-Bodv Velocity Scan Input 



CARD 27: MOLD, MEET, NEAR, INCPRI , INCPRJ, INCPRK 


END OF OFF-BODY VELOCITY 
SCAN INPUT 


MOLD? 
ON CARD 27 


CARD 28: XO, YO, ZO, NP 


Individual YES / NP=1? \ 

' Point \ON CARD 28, 


CARD 29; XI, Yl, Zl, NPl 


NPl < O ?\ 
)N CARD 29, 


CARD 2 9A : (AL(I), 1=1, I NPl| ) 


Points 
Along a 
Line 


YES / NP=2 ? \ 

VN CARD 28. 


Lines 
'within a 
Plane 


CARD 30: X2 , Y2, Z2 , NP2 


NP2<0? \ YES : ^ 

)N CARD 30/ ^ CARD 30A: (AL2(I), 1=1, | NP2 | ) 


YES / NP=3? ' 

\)N CARD 28, 



CARD 31: X3 , Y3, Z3, NP3 


Planes 


/ NP3<0? ' 
ON CARD 3! 


CARD 31A: (AL3(I), 1=1, | NP3 | ) 











ORIGINAL PAGE 18 
OF POOR QUALITY 


OFF-BODY VELOCITY SCAN INPUT, CONTINUED 



12 8 
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S . 0 O UTPUT r>B S CR IPTION 
8 . 1 Prin t Fi le 


The output starts with the program header followed by the 
user's TEXT identification. The list of BASIC INPUT is printed 
next, card by card. Each variable is identified in parentheses. 
The input list for PATCH GEOMETRY is printed next if the IPPI input 
option has been set >0. If IP"I“1 this input list excludes the 
basic point coordinates; these can be included in the list by 
setting IPf»T».2, but this can give a large volume of print for 
complicated configurations. 

If the negative IPPOCM option has been used on CARP 2A, then 
the if ASIC POINTS defining patch SECTIONS are printed patch by 
patch: 


LIST OF BASIC POINTS ON SECTION * OH PATCH fr 

ly dz (he?) 


These basic points include those generated by the code as 
well as those input by the user. The subscripts KBP identify the 
basic points for the purpose of copying strings of points p v e r 
into later sections. 

After processing the patch geometry input, the program 


(subroutine GEOKIN) 
parameters : 

prints 

out a 

s tsB.ary 

table 

of the patch 

Li 

I DENT CLASS 

no; tp 

nnou 

NCOL 

I PAN 

LPAN HP AN 

» 

• e 

• 

• 

• 

» 

• * 

* 

• * 

• 

• 

# 

• 

• e 

• 

* ♦ 

* 

• 

• 

• 

• * 


where I! is the patch number; I PENT is the patch identifier (1, 2 
or 3); IILASS is the ASSEHLLY number and EOKP is the component 
number to which this patch is assigned; NRCV and NCOL are the 
number of panel rows and columns, respectively, on each patch; 
and I PAN, LPAN are the first and last panel subscripts on each 
patch. NPAH is the total number of panels within the patch. If 
a text identification, PNAIIE, has been given for the patch, this 
is written to the right of KPAN, 

If IPEGOK-l (see CARDS 2, under IPKLEV, and/or CARD 2A), 
then the panel corner point coordinates are printed patch by 
patch: 
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PANEL CORNER POINTS OH PATCH # 



Z Z 1 Yi Zj 

* ® • e 

9 ® 9 ® 

« ® 9 ® 


x 2 y 2 z 2 

e e * 

* « 9 

* ♦ • 


j : 3 y 3 z 3 

® • * 

9 ® 9 

• • * 


x 4 ?4 Z 4 

• • • 

* • • 

• • » 


wuerc » k is tne panel subscript# and * j # y^, sj are expressed in 
the global coordinate system. 

If I PECO!! *2 (see CARD 2 under IPELEV and/or CARD 2A) tLen 
tbe panel control point (SC, YC, ZC) and unit normal vector (EM, 
YU, ZU) are printed patch by patch: 


PANEL CONTROL POINTS AMD UNIT NORMAL VECTORS ON PATCE # 


SC 


YC 


ZC 


SN YN EH 


Tl#e list for VAKE INPUT follows. Again, the input variables 
are identified in parentheses. V hen the vfake input Las been 
processed, the wakc-gr id-plane x-stations are printed in 10F1G.4 
format. 

If IPENAr>0 (see CARP 2 under IPELEV and/or GAEL 2A), then the 
panel neighbor information (see Figure 24) is printed patch by 
patch 'before' and 'after wake shedding'. 


PANEL NEIGHBORS ON PATCH # 


LEFCEE UA1IE SEEDLING 
IJADOli (MAESID) 


AFTER VASE SHEETING 
N AEG II (NAB SID) 


EPAN SIDE 1 SIDE 2 SIDE 3 SIDE 4 


EPAN SIDE 1 SIDE 2 SIDE 3 SIDE 4 


If IPUNAE=1 only those panels with missing neighbors (e.g., at 
a wake-shedding line) o$ with special neighbors (e.g., at a plane 
of symmetry) are included in the list. If IFZNAP=2 , all the panel 
neighbor information is printed. 

Each wake shape iteration ( !’.0DE**1 on C ADD 3) starts with the 
printout : 

VASE ITERATION ft 
ALPHA m ft 
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(Unless the print frequency control parameter* IPRX.SS, (see CAR’' 2) 
has suppressed it.) 

On subsequent viscous «- , p otential flow iterations, this v. 7 ill be 
preceded by : 


VISCOUS-'POTEiJTIAL FLOU ITERATION I 

If NODE*2 , the corresponding printout (unless suppressed by 
the parameter IPRESS) is 

IT * * TAU - # 

ALPHA « # 01 .'EGA « # 

where IT is the tine-step number; TAU is the normalized tine from 
the start of the final complete cycle (i.e,, the cycle being 
analysed for real and imaginary terms),* ALPHA is the instantaneous 
incidence (degrees) and ONEGA is the instantaneous rotation rate 
(i.e., ot). 

All the following printout in the wake shape iteration loop 
(liODE=l) or time-step loop (LODE* 8 !) is controlled by the print 
frequency control parameter, IPRESS (CARD 2); the printout can be 
suppressed or just printed for every nth step, say. 

Details of the wake paneling are printed nest depending ox: the 
value of IPRY7AK (see CARD 2 under IPRLEV or CARR 2A); if IPRY.'ADO. 
the following is printed for each w a he: 

LULROUTINL VAN PAN 
LATA FOR T/AEE #, IT = # 

I DEFY? IFLEEV; NY/COL 

* # # 

where IT is the wake shape iteration number (if M C R E = 1 ) or tine 
step nuiaber (if E0DE=2). 

I DEFY.' , IFLEXY'J are the input quantities (see CARD 19) and NY/ COL 
is the number of vale columns on this wake. The text for the wake 
identification (if this option was used ou CARD 19) is printed to 
the right of NYvCCL. 

This is followed by the wake shedding parameters (Figure 27) 
for each column on the wake: 


IDY/COL KYi'PU EY/PUU KUPL EY/PLL SU LSU SL DSL PIUU P II I UUP II ILPE ILL DELV 




where IpVCCL is the wake type (see under CAHIS 19, 2 0); IIVIU, 

Kh'i’iiU, II V PL and KVPLL are subscripts for upper and lower surface 
panels at the wake separation line, see Figure 27; PMIU, HMilMJ, 
P- IL, PiilLL are onset flow velocity potential values at the center 
points of panels KWPU , etc, DELV is the jump in the stro&uwise 
tangential velocity across the wake for those columns with 
IDVCGL=4 ( i . e . , jet wakes). 

If XPr.y AK»2 , the wake line geometry on this wake is also 
printed: 

WAKE LINE GEOMETRY 

(Note: the print excludes the last three points on each line.) 

LINN # 


H X Y 


If IPh.V/.r.=3 , the wake panel doublet values arc printed, column 
by colusn, 

When all the wakes have bees processed, the standard printout 
( i . e . , j£Pl:ttAE-0) from V/ All PAN follows : 

DATA IN ’’.'ANl’AN 

(NV.’C, IV PAN ( IIVC ) , NVC-1, KYCGLT) 

These are printed in sets of 10 across the page. NY.’C is the 
wake column subscript and I VP AN is the wake panel subscript at 
the beginning of each column. NVCOLT is the total number of wake 
columns in the configuration. The wake panel subscripts used in 
the plotting routine CfihIPLT are obtained by adding 3000 to 
IV PAN. 


The relationships between the wake lines and the wake-grid- 
planes is printed next: 

WAKE LIME SECUEIICE ARRAY IN VAIIE-GKID-PLAKE SYSTEM 
J, (LSECUdl, J) II* 1 , NVCP) 


NVGP is the number of w a ke - g r i d-p 1 ane s . J is the wake line 
subscript. 

The w ake-gr id-pl ane s are the columns in the arraj output with 
the upstream plane at the left. Each LCEC.U value is the subscript 
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OF 



of a wake line intersecting each wake-grid-^laue. A negative sign 
is placed on the first and last wake line of each wake . 

A suiiinary of the first wake-grid-plane intersected by each 
wake line is printed next: 

FI -1ST WAKE-GRID-PLANE INTERSECTED BY EACH WAKE LINE 
(J, IV/GP( J ) * J-l, HUT) 
where LINT' is the total number of wake lines. 

Output frot; the doublet solution subroutine LUI-SOL follows 
next. If the cumber of panels is such that the iterative solution 
procedure has been used, then the following is printed. 

First, the blocking arrangement is given, JiFRL is the first 
row of a block, while NROWB is the number of rows in the block. 
This is followed by: 

SOLUTION II I STORY OH PANEL l.'ITIi JlAKII.'UIi RESIDUAL 

JUES1I, SOL (panel subscript and solution value) 
follows in 5(14, 2X, G12.6) format 

If the value of IPRSOL=l (sec under IPRLEV on CARD 2), then 
the complete doublet solution is printed in 5(19* F10.5) format 
(i.e., panel subscript and doublet value). 

The surface pressure distribution is printed next: 

(a) II GDE= 1 . The pressure distribution is printed according to 
the type of patch ( I BELT) . If aDEHT* 1, the following is printed 
for each column of panels. 

If X Y 2 DUE VX VY VZ V C p »C Z^C 


where K is the panel subscript; X, Y, Z is the control point 
location, and DUE the doublet vcluc; VX, VY, VZ and V are the 
velocity vector components and magnitude; C p ^ s t ^ e pressure 
coefficient and X^C, Z^C the location of the po*int relative to the 
local chord line. 


If IPRPPI=1, then the P^P A values are printed to the right of 

Cp- 

Each column is followed by the integrated section force and 
moment coefficients defined in the wind axes. 
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YLr.^DPAi: 


where CD, CS aud CL are the section force coefficients ir. the drag, 
side and lift directions, respectively; ChS, C L Y , C .11 Z the roll, 
pitch and yaw moment coefficients; CHIC is the circulation value; 
CHORD is the local chord length with leading edge at XLE, YLS, ZLE. 
The coordinate are given in the sane reference frame as originally 
specified. YLE-^DSPAK is the normalized spanwise location of the 
section where SSPAH is the reference semispan on CARD 7. SSPA1! is 
used to normalize the rolling and yawing moments while C D A L is 
used to normalize the pitching moment. 


The force and moment coefficients are also given in the body 
axis system. In this case, C II X and C I. Z have been normalized by 
SSPAH and CRY by CBAS. 


For type 2 patches (IDE FT =2), the pressure printout is 
shortened by omitting X-'C, Z^C anu by omitting the section force 
and moment summary at each column. 


A summary of the patch force and moment coefficients 
normalized by the reference quantities in the basic data is 
printed after each patch and the total force and :a o m c n t 
coefficients are printed after the last patch for components and 
assemblies and also for use r- se 1 e c te d panel sets specified under 
the options on CARD 8. 


(b) ?iO DC =2 . The unsteady pressure distribution is printed out 
in a similar way to the RCTZ=1 format. If II.-EuT=i, the following 
is printed for each column of panels on the patch. 


k x y z cpi: cpi cpkcd phase z^c z^c 

• •••• • • • • • 

« » t » • • • * • * 

• •••• • • * • * 


Here, CPR, CPI and CPf.OD are, respectively, the real and imaginary 
pressure coefficients and the pressure coefficient modulus divided 
by ; and PHASE is the phase angle in degrees. Each column is 
followed by a summary of the section real and imaginary lift jand 
pitching moments normalized by local conditions and divided by a. 

CLL CLI CUR CHI 

CL HOD CLPilASE C1.IUOD CMP1IASE 

The phase angles are in degrees. 

This summary and the X-'C , Z*-C details are omitted if IDEKT=2. 

The force and moment summary based on overall conditions is 
printed for each patch and for the total configuration. 
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If I??vCPV = l {see under I PRLEV ou CARD 2 an a*' or CARL- 2 / ) , 
then the pressures a u d velocity results are printed fro a 
subroutine PIIIDIF. 


'Patch No . = . . . ' Patch nutate r of array of corner 

* Column. No, = . . . * point. Col urn u number of array of 

corner points. When an array is 
common to two columns, the snail er 
column number is output. 

NAIM OUTPUT OF CORNER POINT VALUES 


•T> 


z 


PEI VX VY VZ VT CP 


Pill is tie panel doublet value (also the surface perturbation 
potential). IIP is the subscript of the panel containing the 
corner points. When a point is common to two panels the smaller 
panel number is output. L is the location of the corner point 
relative to the panel. Numbering is 1 to 4 , with 1 denoting the 
top left-hand corner. Other corners are numbered sequentially in 
an anticlockwise direction, e.g., see Figure 9. 

After the wake shape iteration loop, the output from 
subroutine STLINE is printed (unless IBLTYP=0 and this is not the 
last viscous -'potential flow solution). 


STLN = n 


*' a ne 1 i < o . — * ... 


bide it o i — ... I* — ... 


u is the stream. line number which is set sequentially by the 
order of data cards input. Other in formation output in the same 
line reproduces the input information on the data card, n, except 
for Side No. If NS was set to zero in the input, then Side No. 
is the value determined by the program. 

The main output of points on the streamline are: 

up x y z v:: vy vz vt cp ds ge 


where HP is the panel number; X, Y, Z the point coordinates; VX, 
VY, VZ, VT is the velocity vector and its magnitude; CP is the 
pressure coefficient; DS is the distance of the point along the 
streamline; and CX has the significance of geodesic curvature on 
a body. Cm a wing it serves as a measure of convergence or 
divergence of s tr earal ine s . 

The output from the boundary layer routine follows if NVPI>0 
ou CARD 4(a). The calculated boundary layer characteristics arc 
tabulated for each streamline " and include the local shape 
parameter and shin friction coefficient. Each table is preceded 
by a summary of the streamline geometry and pressure 
distribution. 


136 



Output f roc tke off-body velocity scan routine coses nest. 
iOr each scan box the input data is printed first, followed by 

I Ml' OX MOUNT 


where I is the scan-plane subscript in the box, MDOX; HOUIIT is a 
number written to the plot file (see below) to identify the scar, 
plane. The velocity and pressure data is then printed point by 
point for each line in the plane. 


VX 


VY 


VZ 


C£ 


This data is also written to tke plot file (TAPE 7) except that 
XCIiMT is substituted for X. 


Output from the off-body streamline routine comes last. The 
input data for each streamline is followed by the streamline 
geometry and aerodynamic data. 


VX 


VY 


VZ 


Co 


here, I! is the local Mack number and S is the distance 
from the upstream end. 


ensured 


8.2 Plpt_File 

During the execution of the program, an unformatted plot 
file (TAPE 7) is assembled as follows: 


(i) Our face Panel Geometry 

For each panel there are four records containing the 
four corner points. 


UP, 

V «| 

Yl, 

Z1 

NP. 

12 , 

Y2 , 

Z2 

UP, 

S3, 

Y3, 

Z3 

UP, 

X4, 

Y4 , 

Z4 


There are Ms sets. 
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ORIGINAL PAGE IS 
OF POOR QUALiTf 

Cii) y a 1: e Panel Geometry 

The v/&ke panel corner points are written in the same ray as 
the surface panels except the counter IIP starts at 3001. A f ter 
each pass the wake panel geometry is terminated by 

-IT, 0.0, 0.0, 0.0 

where IT is the wake shape iteration number (!!0DE*1) or the time- 
step number (I101>£*2) . 

For each cycle the aerodynamic data follows this line: 

Nf, (the total number of panels) 

followed by Ns records: 

HP, ZC, YC, ZC, VS, VY, VZ, V, CP 

If I10PE = 2, the last set of aerodynamic data is followed by 
the unsteady pressure analysis; i.e., N s records of 

HP, XC, YC, ZC, CPr, CPl, CPllOD* PllASE 


Ofi’-body velocity data is placed next on the file 

i:ouNT, x, y. z, vx* vy, vz, v, cr 

where ECUKT starts at <5001 and is incremented by 1 for each scan 
line, 1C for a scan plane and ICC for a box. 

- 1 , 0 . 0 , 0 . 0 , 0.0 

The streamline data comes next. The first record is 

KLIIIES 

This is the total number of streamlines (i.e., on-body and off- 
body). This is followed by one record for each streamline. 

I LAX , ( (ST( J,E) , 1*1, 12), KPRL(X) , X=L, IIIAX) 

where the quantities in the ST array include 

X, Y, Z, VX, VY, VZ, V, CP. S. 

The HPSL(Z) values are surface panel subscripts crossed by the 
streamline. Off-body streamlines are treated in a similar way 
bat HP SL values are zero in this case. The streamline data is 
terminated with a -1 for HPSL. 
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Finally, the boundary layer data is placed on the tape with 
two records per surface streamline, 

JLINE - the streamline number ; and 

KPOINT , (D (f ) , 1=1, 13) 

The records are terminated by setting JLINE®-!. 

The boundary layer parameters are: 

D(l) * x-coordinate 
0(2) ■ y-coordina te 
0(3) » z-coordinate 

0(4) ■ s , distance along streamline 
D(5) * v, velocity 

D(6) “ 9u^3s , rate of change of velocity along s direction 
D(7) * H , shape factor 
D(8) ■ 6, displacement thickness 
D(9) ■ THT, momentum thickness 
D(10) * RTH, Reynolds number based on momentum thickness 
0(11) * CFD , skin friction drag 

D(12) = CDS, skin friction from Squire and Young formula 
D(13) * Displacement thickness divided by y-coordinate 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


9.0 

1 • M» skew , B . , 1 Pr e diction of Subsonic Aerodynamic Character- 
istics: A Case for Low-Order Panel Methods ' , £. Aircraft . 

Vol . 19, No. 2 , Febrnary 1982 , pp . 157-163 . 

2 . Maskew, B . , Rao, B.M and Dvorak, F.A., 'Prediction of Aero- 
dynamic Characteristics for lings with Extensive Separations', 
Paper No. 31 in Connotation &f Yi#cons-Inyis c i d Interactions, 
AGARD-CP-291 , February 1981. 

3. Maskew, B. , 'Influence of Rotor Blade Tip Shape on Tip Vortex 
Shedding — An Unsteady Inviscid Analysis', Proceedings of the 
36th Annual Forum of the American Helicopter Society, Paper 
80-6, May 1980. 

4. Sanaa, J.M. and Maskew. B. , 'A Surface Singularity Method for 
Rotors in Hover or Climb', USAAVRADCOM-TR-81-D-23 ■ December 
1981. 

5. Clark, D.R. and Maskew, B . , 'An Analysis of Airf rame^Rotor 
Interference in Forward Flight', Paper No. 50, Presented at 
the 7th European Rotor craft and Powered Lift Aircraft Forum, 
Garmisch-Partenkirchen, FRG, September 1981. 

6. Maskew, B . , 'A Quadrilateral Vortex Method Applied to Con- 
figurations with High Circulation', Paper No. 10, Presented 
at the NASA Workshop on Vortex-Lattice Utilization, Langley 
Research Center, May 1978, Also NASA SP-405 , May 1980. 

7. Dvorak, F.A. and Woodward, F . A . , ' Vi scou s ^Potent ial Flow An- 

alysis Method for Multi-Element Infinite Swept Wings', NASA 
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8. Dvorak, F.A., Maskew, B. and Woodward, F.A. , 'Investigation of 
Three-Dimensional Flow Separation on Fuselage Configurations', 
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Research and Development Laboratory, Fort Eustis, VA 23604, 
March 1977. 
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J.i tiie following gages the input and output art „iven tor the 
\i i a g - 1) o u 7 totf ijjuration sLovu in Figure AI . Tlie way the surlace 
ha - been broken into seven patches is shown it: Figure A 2 (a). Each 
patch is iliustratec separately i a Figures k 2 { b > t h ? o u „ h { 1: ) , 
si. cm i n A the form of input. 

The figures are followed by a listing of the coEflite i up ut 
etch. Finally, samples frou the program output file arc given. 


14 
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Figure ftl. Wing-Body Sample Case 
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utxc-iODv AC»oO'VM*niCS,TE$T CASt Figure A2(b). Patch 1, Wing. 





Q 
*— * 
CO 
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SEMI -ELLIPSE 



u!NG-tODV AERODVNMHCS'TCST CftSC Figure A2 (d) . Patch 3, Forebody 










NOSEC * 34 



14? 









NDSEC = 42 (COPIED FROM WING) 
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INPUT LISTING 


WING/BODY TEST CASE 
1 0 1 

1 SCO 
0 1 1 


0 J 

1C 0 


£ 05 

145 2 

12 0 

2 

304 

3 26 

1 

331 

3 1 

1 

1 

0 1 

1 WING 

14 325 

1. 6 

0. 0 

100. 0 

0 0 

. 8592 

85 0 

- 865 

. 9718 

70 0 

-1. 935 

. 9718 

55. 0 

-2 74 1 

. 9310 

40 0 

-2. 998 

. 662 

30. 0 

-2. 852 

. 507 

20 0 

-2. 482 

. 3099 

10 0 

-1 . 824 

. 0634 

5 0 

-1 310 

0704 

2 5 

-. 956 

-. 1408 

1 25 

-. 717 

169 

0 ^5 

574 

-. 1972 

0 5 

-. 476 

— , 2183 

0 0 

0. 0 

~. 2254 

CL 

476 

-. 2183 

0. 75 

574 

1972 

1 . 25 

717 

-. 169 

2 5 

. 956 

- 1 4C8 

5. 0 

1. 310 

0704 

10 0 

1. 824 

. 0634 

20. 0 

2. 482 

. 3099 

30. 0 

2 852 

. 507 

40. 0 

2. 998 

. 662 

55 0 

2. 741 

. 9310 

70 0 

1 935 

.9718 

es o 

0. B65 

. 9718 

100. 0 

0. 0 

. 8592 

25. 375 

12 0 

0. 0 

100. 0 

0. 0 


85 0 

-. 865 


i 70. o 

- 1 . 935 


I 55. 0 

-2. 741 


j 40. 0 

-2. 998 


! 30. 0 

-2 852 


j 20. 0 

-2 482 


| 10 0 

-1. 824 


| 5. 0 

-1 310 


( C. 

C_ . w' 

-. 956 


1. 25 

- 717 


75 

574 


u. 5 

476 


| 0 0 

J 

0 0 


0 5 

476 


j 0 75 

. 574 



0 0 


. 071 


1 


3 

. 045 


1 


0. 0 


0. 0 





0. 0 


O. 0 


2 


3 5 



: £5 

. 717 

2. 5 

. 956 

5 0 

1. 310 

10 0 

1. 824 

0 

2. 482 

30. 0 

2. 852 

40. 0 

2. 998 

*5 0 

2. 741 

70. 0 

1 935 

85, 0 

0. 865 

100 0 

0. 0 



c. 

1 

1 

1 WING 

TIP 

3 3 

0 


2 

0 

1 

1 FOREBODY 




1 o 
»-* ' 
o 
o 

O 

6 

1 

0 0 

Q 

6 

o 

6 

o 

6 

6 

o 

o 


6 


3 




2 0 

0. 0 


0. 0 

. 4830 

0, 0 

0. 0 

0 

4 0 

0. 0 


0. 0 

. 7900 

0. 0 

0. 0 

0 

6. 0 

0. 0 


0. 0 

1. 0775 

0. 0 

0. 0 

0 

8 0 

. 0 
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16 Abstract 

VSAER0 is a computer program used to predict the non-linear aerodynamic 
characteristics of arbitrary three dimensional configurations in subsonic 
flow. Non-linear effects of Vortex Separation and Vortex/ Surface inter- 
action are treated in an iterative wake-shape calculation procedure, while 
the effects of viscosity are treated in an iterative loop coupling 
potential-flow and integral boundary-layer calculations. The program 
employs a surface singularity panel method using quadrilateral panels on 
which doublet and source singularities are distributed in a piecewise 
constant form. This User's Manual provides a brief overview of the mathe- 
matical model, instructions for configuration modelling and a description 
of the input and output data. A listing of a sample case is included. 
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